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Abstract 



After the separation of the center-of-mass motion, a new privileged class of 
canonical Darboux bases is proposed for the non-relativistic N-body problem 
by exploiting a geometrical and group theoretical approach to the definition of 
body frame for deformable bodies. This basis is adapted to the rotation group 
S0(3), whose canonical realization is associated with a symmetry Hamiltonian 
left action. The analysis of the S0(3) coadjoint orbits contained in the N-body 
phase space implies the existence of a spin frame for the N-body system. Then, 
the existence of appropriate non-symmetry Hamiltonian right actions for non- 
rigid systems leads to the construction of a N-dependent discrete number of 
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dynamical body frames for the N-body system, hence to the associated notions 
of dynamical and measurable orientation and shape variables, angular velocity, 
rotational and vibrational configurations. For N=3 the dynamical body frame 
turns out to be unique and our approach reproduces the xxzz gauge of the 
gauge theory associated with the orientation- shape SO (3) principal bundle 
approach of Littlcjohn and Rcinsch. For > 4 our description is different, 
since the dynamical body frames turn out to be momentum dependent. The 
resulting Darboux bases for A'^ > 4 are connected to the coupling of the spins 
of particle clusters rather than the coupling of the centers of mass (based 
on Jacobi relative normal coordinates). One of the advantages of the spin 
coupling is that, unlike the center-of-mass coupling, it admits a relativistic 
generalization. 
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I. INTRODUCTION. 



This paper deals with the construction of a speciahzed system of coordinates for the non- 
relativistic N-body problem which could be instrumental to nuclear, atomic and molecular 
physics, as well as to celestial mechanics. In particular, we shall exploit the technique of 
the canonical realizations of Lie symmetry groups within the framework of the non- 

relativistic version of the rest-frame Wigner covariant instant form of dynamics to the 
effect of obtaining coordinates adapted (locally in general) to the SO (3) group. In most of 
the paper we consider only free particles, since the mutual interactions are irrelevant to the 
definition of the kinematics in the non-relativistic case. 

Isolated systems of N particles possess 3N degrees of freedom in configuration space 
and 6N in phase space. The Abelian nature of the overall translational invariance, with 
its associated three commuting Noether constants of the motion, allows for the decoupling 
and, therefore, for the elimination of either three configurational variables or three pairs of 
canonical variables, respectively (separation of the center-of-mass motion). In this way one 
is left with either 3N-3 relative coordinates pa or 6N-6 relative phase space variables p^, p^, 
a = 1, .., — 1 and the center-of-mass angular momentum or spin is = Y^a=i Pa x Pa- In 
the non-relativistic theory most of the calculations employ the sets of 3A^ — 3 Jacobi normal 
relative coordinates Sa which diagonalize the quadratic form associated with the relative 
kinetic energy [the spin becomes S = Z^^i^ Sa x vf^a, with vf^a momenta conjugated to the 
Sas\. Each set of relative Jacobi normal coordinates Sa, a = 1, .., N — 1, is associated with 
a different clustering of the N particles, corresponding to the centers of mass of the various 
subclusters. In special relativity Jacobi normal coordinates do not exist, as it will be shown 
in Ref. 0, and a different strategy must be used. 

On the other hand, the non-Abelian nature of the overall rotational invariance entails 
the impossibility of an analogous intrinsic separation of rotational (or orientational) config- 
urational variables from others which could be called shape or vibrational. As a matter of 
fact, this is one of the main concerns of molecular physics and of advanced mechanics of 
deformable bodies. Recently, a new approach inspired by the geometrical techniques of fiber 
bundles has been proposed in these fields of research: a self-contained and comprehensive 
exposition of this viewpoint and a rich bibliography can be found in Ref. 0. 

In the theory of deformable bodies one looses any intrinsic notion of body frame, which is 
a fundamental tool for the description of rigid bodies and their associated Euler equations. 
A priori, for a given configuration of a non-relativistic continuous body, and in particular 
for a N-body system, any barycentric orthogonal frame could be named body frame of the 
system with its associated notion of vibrations. 

This state of affairs suggested to replace the kinematically accessible region of the non- 
singular configurations Q in the (3N-3)-dimensional relative configuration space by a SO (3) 



^ See Refs. [p|,[lO| for a discussion of the singular (collinear and N-body collision) configurations. 
Applying the SO (3) operations to any given configuration of the 3N-3 relative variables (a point 
in the relative configuration space) gives rise to 3 possibilities only: i) for generic configurations 
the orbit containing all the rotated copies of the configuration is a 3-dimensional manifold [diffeo- 
morphic to the group manifold of S0(3)]; ii) for collinear configurations the orbit is diffeomorphic 



3 



principal fiber bundle over a (3N-6)-dimensional base manifold, called shape space 0. The 
S0(3) fiber on each shape configuration carries the orientational variables (e.g. the usual 
Euler angles) referred to the chosen body frame. A local cross section of the principal fiber 
bundle selects just one orientation of a generic N-body configuration in each fiber (SO (3) 
orbit) and this is equivalent to a gauge convention., namely to a possible definition of a body 
frame {reference orientation)., to be adopted after a preliminary choice of the shape variables. 
It turns out that this principal bundle is trivial only for N=3, so that in this case only global 
cross sections exist, and in particular the identity cross section may be identified with the 
space frame. In this case any global cross section is a copy of the 3-body shape space and 
its coordinatization gives a description of the internal vibrational motions associated with 
the chosen gauge convention for the reference orientation. For > 4, however, global cross 
sections do not exist and the definition of the reference orientation {body frame) can be 
given only locally. This means that the shape space cannot be identified with a (3N-6)- 
dimensional submanifold of the (3N-3)-dimensional relative configuration space. The gauge 
convention about the reference orientation and the consequent individuation of the internal 
vibrational degrees of freedom requires the choice of a connection on the SO (3) principal 
bundle (i.e. a concept of horizontality) and this leads in turn to the introduction of a SO (3) 
gauge potential on the base manifold. Obviously, physical quantities like the rotational or 
vibrational kinetic energies and, in general, any observable feature of the system must be 
gauge invariant, namely independent of the chosen convention. Note that both the space 
frame and the body frame components of the angular velocity are gauge quantities in the 
orientation- shape bundle approach and their definition depends upon the gauge convention. 

While a natural gauge invariant concept of purely rotational N-body configurations exists 
(when the N-body velocity vector field is vertical, i.e. when the shape velocities vanish), a 
notion of horizontal or purely vibrational configuration requires the introduction a connection 
r on the SO (3) principal bundle. A gauge fixing is needed in addition in order to select a 
particular F-horizontal cross section and the correlated gauge potential on the shape space. 
See Ref. for a review of the gauge fixings used in molecular physics' literature and, in 
particular, for the virtues of a special connection C corresponding to the shape configurations 
with vanishing center-of-mass angular momentum S 0. 

This orientation- shape approach replaces the usual Euler kinematics of rigid bodies and 
entails in general a coupling between the internal shape variables and some of the orienta- 
tional degrees of freedom. In Ref. it is interestingly shown that the non-triviality of the 



to tlie 2-sphere 5^; iii) for tfie N-body collision configuration (in whidi all tfie particles coincide at 
a single point in space) the orbit is a point. 

^It is the open set of all the orbits of generic non-singular configurations. 

^This is due to the topological complexity of the shape space generated by the singular configu- 
rations which are dispersed among the generic configurations for > 4. 

^The C-horizontal cross sections are orthogonal to the fibers with respect to the Riemannian 
metric dictated by the kinetic energy. 
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so (3) principal bundle, when extended to continuous deformable bodies, is at the heart of 
the explanation of problems like the falling cat and the diver. A characteristic role of SO (3) 
gauge potentials in this case is to generate rotations by changing the shape. 

In Ref. [§] the Hamiltonian formulation of this framework is also given, but no explicit 
procedure for the construction of a canonical Darboux basis for the orientational and shape 
variables is worked out. See Refs. for the existing sets of shape variables for = 3,4 

and for the determination of their physical domain. 

Independently of this SO (3) principal bundle framework and having in mind the rel- 
ativistic N-body problem where only Hamiltonian methods are available, we have been 
induced to search for a constructive procedure for building canonical Darboux bases in the 
(6N-6)-dimensional relative phase space, suited to the non-Abelian canonical reduction of 
the overall rotational symmetry. Our procedure surfaced from the following independent 
pieces of information: 

A) In recent years a systematic study of relativistic kinematics of the N-body problem, 
in the framework of the rest-frame Wigner covariant instant form of dynamics has been 
developed in Ref. [|r^ and then applied to the isolated system composed by N scalar charged 
particles plus the electromagnetic field |[T^JT3[| . 

These papers contain the construction of a special class of canonical transformations, 
of the Shanmugadhasan type [1^J§]. These transformations are simultaneously adapted to: 
i) the Dirac first class constraints appearing in the Hamiltonian formulation of relativistic 
models; ii) the timelike Poincare orbits associated with most of their configurations. In the 
Darboux bases one of the final canonical variables is the square root of the Poincare invariant 

(P^ is the conserved timelike four- momentum of the isolated system). Subsequently, by 
using the constructive theory of the canonical realizations of Lie groups a new family 
of canonical transformations was introduced in Ref. IITHI. This latter leads to the definition 



of the so-called canonical spin bases, in which also the Pauli-Lubanski Poincare invariant 
W'^ = —P'^S^ for timelike Poincare orbits becomes one of the final canonical variables. 
The construction of the spin bases exploits the clustering of spins rather than the Jacobi 
clustering of centers of mass. 

In spite of its genesis in a relativistic context, the technique used in the determination of 
the spin bases, related to a typical form of the canonical realizations of the E(3) group, 
can be easily adapted to the non-relativistic case, where W"^ is replaced by the invariant S"^ 
of the extended Galilei group. 

B) These results provide the starting point for the construction of a canonical Darboux 
basis adapted to the non-Abelian SO (3) symmetry. The three non-Abelian Noether con- 
stants of motion 5* are arranged in these canonical Darboux bases as an array containing 
the canonical pair S^, P = tg~^^ and the unpaired variable S = \S\^ (scheme A of the 
canonical realization of SO (3) [Q). The angle canonically conjugated to S, say a, is an 
orientational variable, which, being coupled to the internal shape degrees of freedom, cannot 



^For the configurations of the isolated system having a rest-frame Thomas canonical spin St 
different from zero. 

^In this context the configurations with S = are singular and have to be treated separately. 
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be a constant of motion. In conclusion, in this non-Abelian case one has only two (instead 
of three as in the Abelian case) commuting constants of motion, namely S and (like in 
quantum mechanics). 

This is also the outcome of the momentum map canonical reduction [p!^ -[l8| by means of 
adapted coordinates. Let us stress that a, S^, [3 are a local coordinatization of any coadjoint 
orbit of SO (3) contained in the N-body phase space. Each coadjoint orbit is a 3-dimensional 
embedded submanifold and is endowed with a Poisson structure whose neutral element is 
a. This latter is also the essential coordinate for the definition of the flag of spinors (see 
Refs. I TP] and [Q, Section V). On the other hand, the spinor flag is nothing else than a unit 
vector R orthogonal to S [^, which is going to be a fundamental tool in what follows. 

By fixing non-zero values of the variables S''^, (3 = tg~^^ through second class constraints, 
one can define a (6N-8)-dimensional reduced phase space. However, the canonical reduction 
cannot be furthered by eliminating S, since a is not a constant of motion. 

C) The group-theoretical treatment of rigid bodies [Chapter IV, Section 10] is based 
on the existence of the realization of the (free and transitive) left and right Hamiltonian 
actions of the SO (3) rotation group on either the tangent or cotangent bundle over their 
configuration space. Given a laboratory or space frame fr, the generators of the left Hamilto- 
nian actio nQ are the non-Abelian constants of motion S^, S"^, S^, [{S''', S'^} = e''"'*'"S'"], viz. the 
spin components in the space frame. In the approach of Ref. the SO (3) principal bundle 
is built starting from the relative configuration space and, upon the choice of a body-frame 
convention, a gauge-dependent SO (3) right action is introduced. 

Similarly, taking into account the relative phase space of any isolated system, one may 
investigate whether one or more SO (3) right Hamiltonian actions could be implemented 
besides the global canonical realization of the SO (3) left Hamiltonian action, which is a 
symmetry action. In other words, one may look for solutions 5"'", r=l,2,3, [with J^riS^Y = 
^^(5'-)2 = S^], of the partial differential equations {S',S'} = 0, {S',S'} = -e^'^"5" and 
then build corresponding left invariant Hamiltonian vector fields. Alternatively, one may look 
for the existence of a pair 5*^, 7 = tg~^^, of canonical variables satisfying {7,5*^} = — 1, 
{7, S'^} = {S^, S^} = and also {7, a} = {S^, a} = 0. Local theorems given in Refs. [|l],|^ 
guarantee that this is always possible provided > 3. See Chapter IV of Ref. for 
what is known in general about the actions of Lie groups on symplectic manifolds. Clearly, 
the functions do not generate symmetry actions, because they are not constants of the 
motion. 

The inputs coming from A), B), C) together with the technique of the spin bases of Ref. 
||T5| suggest the following strategy for the geometrical and group-theoretical identification 
of a privileged class of canonical Darboux bases for the N-body problem: 

1) Every such basis must be a scheme B (i.e. a canonical completion of scheme A) 
for the canonical realization of the rotation group SO (3), viz. it must contain its invariant 
5* and the canonical pair 5^, (3 = tg~^^. Therefore, all the remaining variables in the 



^We follow the convention of Ref. |^] ; note that this action is usually denoted as a right action in 
mathematical texts. The S^'s are the Hamiltonians, associated with the momentum map from the 
symplectic manifold to so(3)* [the dual of the Lie algebra so(3)], which allow the implementation 
of the symplectic action through Hamiltonian vector fields. 
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canonical basis except a are SO (3) scalars. 

2) As said above, the existence of the angle a satisfying {a, S*} = 1 and {a, S^^ = 
{a, (3} = leads to the geometrical identification of a unit vector R orthogonal to S and, 
therefore, of a orthonormal frame S, R, S x R, which will be called spin fram^. 

3) The study of the equations R^ = 1 and {S ■ R, R^} = entails the symplectic result 
{i?*, R^} = 0. As a byproduct we get a canonical realization of an E(3) group with generators 
S, R [{R\ W} = 0, {R\ = e^i^R^] and fixed values of its invariants R^ = l,S-R = Q{ 
non- irreducible type 3 realization according to Ref. [p!5[] ). 

4) In order to implement a SO (3) Hamiltonian right action in analogy with the rigid 
body theory |T^, we must construct an orthonormal triad or body frame N, x, N x x- The 



decomposition 

S = S^X + S^N XX + S^N =^ S''er, (1.1) 

identifies the SO (3) scalar generators S**" of the right action provided they satisfy {S''^, S''^} = 
_^rsugu^ This latter condition together with the obvious requirement that iV, x, iV x x be 
S0(3) vectors [{iV^ 5^} = e^'^^iV", {x\S'} = e"^^", {N x x , S'} = e""iV x £] entails 
the equations [| 

{iv^ N'} = {iv^ x'} = {x^ r } = o. (1.2) 

To each solution of these equations is associated a couple of canonical realizations of the 
E(3) group (type 2, non- irreducible): one with generators S, N and non- fixed invariants 

= S ■ N and |A^|; another with generators S, x ^md non-fixed invariants = S ■ x 
and Ixl- These latter contain the relevant information for constructing the angle a and 
the new canonical pair S^, 7 = ^5'~^fT of S0(3) scalars. Since {a,S^} = {0,7} = must 
hold, it follows that the vector A^ necessarily belongs to the S-R plane. The three canonical 
pairs S, a, S^, (3, S^, 7 will describe the orientational variables of our Darboux basis, while 
|A^| and \x\ will belong to the shape variables. Alternatively, an anholonomic basis can be 
constructed by replacing the above six variables by S**" and three uniquely determined Euler 
angles a, /3, 7. 

Let us consider the case N=3 as a first example. It turns out that a solution of Eqs.( [1.2|) 
corresponding to a body frame determined by the 3-body system configuration only, as in 
the rigid body case, is completely individuated once two orthonormal vectors A^ and Xi 
functions of the relative coordinates and independent of the momenta, are found such that 
A^ lies in the S-R planeQ. We do not known whether in the case N=3 other solutions 



°The notation ' means unit vector. 

^With S'' = S-er, the conditions {5^ S'} = -e"''^" imply the equations S-Cr x Cs+S^S^el, 4} = 
^rsuS^e"^, hence the quoted result. 

^'^Let us remark that any pair of orthonormal vectors N, x function only of the relative coordinates 
can be used to build a body frame. This freedom is connected to the possibility of redefining a body 
frame by using a configuration-dependent arbitrary rotation, which leaves in the S-R plane. 
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of Eqs. ( |1.2|) exist leading to momentum dependent body frames. Anyway, our constructive 
method necessarily leads to momentum-dependent solutions of for Eqs. (|1.2| ) for > 4 and 
therefore to momentum-dependent or dynamical body frames. 

We can conclude that in the N-body problem there are hidden structures allowing the 
identification of special dynamical body frames which, being independent of gauge conditions, 
are endowed with a physical meaning. 

The following particular results can be proven: 

i) For = 2, a single E(3) group can be defined: it allows the construction of an 
orthonormal spin frame S , R, R x S in terms of the measurable relative coordinates and 
momenta of the particles. 

ii) For N = 3, S = Si + S2, a. pair of E(3) groups emerge associated with Si and 5^2, 
respectively. In this case, besides the orthonormal spin frame, an orthonormal dynamical 
body frame N, x, N x x can be defined such that S^ = S ■ x, S'^ = S ■ N x x, S^ = S ■ N are 
the canonical generators of a SO (3) Hamiltonian right action. The non-conservation of S^ 
entails that the dynamical body frame evolves in a way dictated by the equations of motion, 
just as it happens in the rigid body case. 

It will be shown that for N=3 this definition of dynamical body frame can be reinterpreted 
as a special global cross section (xxzz gauge, where x stays for x and z for A^; this outcome 
is independent from the particular choice made for A^ and x) of the trivial SO (3) principal 
bundle of Ref. namely a privileged choice of body frame. Actually, the three canonical 
pairs of orientational variables S^, (3 = tg~^^; S, a; S^, 7 = tg~^^, can be replaced by 

the anholonomic basis of three Euler angles a, (3, 7 and by S^, S"^, S^ as it is done in Ref. 
P]. In our construction, however, the Euler angles a, (3, 7 are determined as the unique set 
of dynamical orientation variables. Then, the remaining canonical pairs q^, p^, /i = 1,2,3, 
of this spin-adapted Darboux basis describe the dynamical shape phase space. 

While the above dynamical body frame can be identified with the global cross section 
corresponding to the xxzz gauge, all other global cross sections cannot be interpreted as 
dynamical body frames (or dynamical right actions), because the SO (3) principal bundle of 
Ref. is built starting from the relative configuration space and, therefore, it is a static, 
velocity-independent, construction. As a matter of fact, after the choice of the shape con- 
figuration variables q^ and of a space frame in which the relative variables have components 
p^, the approach of Ref. begins with the definition of the body-frame components Pa{q^) 
of the relative coordinates, in the form pj^ = R^^{9°')pl^{q^) |^, and then extends it in a 

def 

velocity-independent way to the relative velocities ff^ = i?''*(6'")t'^. In our construction we 
get instead p"^ = R^^{a, j3,^)pl^{q^) in the xxzz gauge, so that in the present case (N=3) all 
dynamical variables of our construction coincide with the static variables in the xxzz gauge. 
On the other hand, in the relative phase space, the construction of the evolving dynamical 
body frame is based on non-point canonical transformations 

iii) For N=4, where S = Si + S2 + S3, it is possible to construct three sets of spin 
frames and dynamical body frames corresponding to the hierarchy of clusterings {{ab)c) 



^^R is a rotation matrix, 0" are arbitrary gauge orientational parameters and p\{q^) is assumed 
to depend on the shape variables only and not on their conjugate momenta. 
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[i.e. ((12)3), ((23)1), ((31)2)] of the relative spins Sa The associated three canonical 
Darboux bases share the three variables S^, (3, S (viz. S), while both the remaining three 
orientational variables and the shape variables depend on the spin clustering. This entails 
the existence of three different SO (3) right actions with non-conserved canonical generators 
Sl^A)y A=l,2,3. Therefore, one can define three anholonomic bases a{A), P{A), 7(A), S^a) 
associated shape variables q^^y P(A)fj., = l,--,6, connected by canonical transformations 
leaving fixed. These anholonomic bases and the associated evolving dynamical body 
frames, however, have no relations with the N=4 static non-trivial SO (3) principal bundle 
of Ref. [§[, which admits only local cross sections. As a matter of fact, one gets = 

R'^'iaiA), f3{A),l{A))piA)a{q(A)^P{A),i, SIa)) mstcad of = i?"(a, ^, 7)/5^(g^). 

These results imply that, for N=4, the 18-dimensional relative phase space admits three 
operationally well defined dynamical body frames, and associated right actions, and its co- 
ordinates are naturally splitted in three different ways into 6 dynamical rotational variables 
and 12 generalized dynamical shape variables. As a consequence, we get three possible 
definitions of dynamical vibrations. Each set of 12 generalized dynamical canonical shape 
variables is obviously defined modulo canonical transformations so that it should even be 
possible to find local canonical bases corresponding to the local cross sections of the N=4 
static non-trivial SO (3) principal bundle of Ref. [Q. 

Our results can be extended to arbitrary N, with S = Y.a=i ^a- There are as many 
independent ways (say K) of spin clustering as in quantum mechanics. For instance for 
N=5, K = 15 : 12 spin clusterings correspond to the pattern {{{ab)c)d) and 3 to the pattern 
{{ab){cd)) [a,b,c,d = 1,..,4]. For N=6, K = 105: 60 spin clusterings correspond to the 
pattern {{{{ab)c)d)e), 15 to the pattern {{{ab){cd)e) and 30 to the pattern {{{ab)c){de)) 
[a,b,c,d,e = 1,..,5]. Each spin clustering is associated to: i) a related spin frame; ii) a 

related dynamical body frame; iii) a related Darboux spin canonical basis with orientational 

52 

variables 5*^, (3, S, q;(a), ^^a)^ 7(A) = tg^^ jr^, ^ = 1, [their anholonomic counterparts 

are a(A), /3(a), 7(A), S^a) with uniquely determined orientation angles] and shape variables 
9(A)' Pm(A), a* = 1, ••, 3iV — 6. Furthermore, for > 4 we find the following relation between 
spin and angular velocity: S"' = X''''(g['^)) uf^^ + f''{q'(A))P{A)fi- 
Let us conclude this Introduction with some remarks. 

The S = 0, C-horizontal, cross section of the static S0(3) principal bundle corresponds 
to N-body configurations that cannot be included in the previous Hamiltonian construction 
based on the canonical realizations of SO (3): these configurations (which include the singular 
ones) have to be analyzed independently since they are related to the exceptional orbit of 
SO (3), whose little group is the whole group. 

While physical observables have to be obviously independent of the gauge-dependent 
static body frames, they do depend on the dynamical body frame, whose axes are operationally 
defined in terms of the relative coordinates and momenta of the particles. In particular, a 



In a way analoguous to the angular momentum composition in quantum mechanics; note that 
this spin clustering is independent of the center-of-mass clustering associated with the Jacobi 
coordinates: the existence of these two unrelated clusterings might prove to be a useful and flexible 
tool in molecular physics. 
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dynamical definition of vibration, wliicli replaces tlie S = C-horizontal cross section of the 
static approach^, is based on tlie requirement tfiat the components of the angular velocity 
vanish. Actually, the angular velocities with respect to the dynamical body frames become 
now measurable quantities, in agreement with the phenomenology of extended deformable 
bodies (e.g. the treatment of spinning stars in astrophysics). 

In Section II the rest-frame description of N free particles together with the Jacobi 
normal relative coordinates is introduced and some further informations are summarized 
about the orientation-shape SO (3) principal bundle approach, both from the Lagrangian 
and the Hamiltonian point of view. 

In Section III the canonical spin bases, the spin frame and the dynamical body frames 
are introduced and the cases = 2, = 3 and A^ > 4 are analyzed separately. 

In Section IV a short account is given of N particles interacting through a potential. 

Some final remarks are given in the Conclusions. 

Appendix A contains the Lagrangian and Hamiltonian equations of motion in the static 
orientation- shape bundle approach. 

In Appendix B the Lagrangian and Hamiltonian results of Subsection D of Section II are 
reformulated in arbitrary (but not Jacobi normal) relative coordinates. 

In Appendix C detailed calculations are given for the case N=3. 

In Appendix D some notions on Euler angles are reviewed. 

In Appendix E the gauge potential in the xxzz gauge is evaluated. 

In Appendix F the construction of the canonical spin bases is given for the A^ = 4 case. 



^Being connected to the Riemannian metric of the non-relativistic Lagrangian, this concept does 
not survive the transition to special relativity anyway. 
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II. THE CENTER OF MASS, THE JACOBI RELATIVE COORDINATES AND 

THE ORIENTATION-SHAPE BUNDLE. 



In this Section we first formulate tlie description of N non-relativistic free particles in the 
rest frame and then we review the theory of the orient at ion- shape principal SO (3) bundle 
of Ref. §. 

Given the coordinates fji{t), i = 1,..,N, of N particles of mass mj, the standard La- 
grangian is L = Y^iLi ^Vi ■ By introducing the center-of-mass coordinates Qnr and a set of 
relative variables, the Lagrangian can be rewritten as L = Y^nr~^ (quadratic form in the 
relative velocities), M = Y^^=i^i- The canonical momenta are pi = miff^, while the total 

momentum conjugated to qnr is P = Y^iLiPi ■ The Hamiltonian is H = YliLi ^ = IM'^ 
(quadratic form in the relative momenta). 

The rest-frame description (equivalent to the decoupling of the center of mass) is obtained 
by imposing the vanishing of the conserved total momentum P = X^ili Pi = 



A. The non-relativistic rest-frame description. 

The rest-frame description of the relative motions can be obtained as the non-relativistic 
c — > cxD limit of the relativistic rest-frame instant form of Ref. |^ . Equivalently we can start 
from the Lagrangian 

Loit) = E -^im + A(t)]^ Sn = J dtLnit). (2.1) 



i=l 



in which the Lagrange multipliers X{t) are considered as configurational variables. 
The canonical momenta are 

Ut) = ^ = 0. (2.2) 

dX{t) 

Therefore, 7rx{t) ^ is a primary constraint. The canonical and Dirac Hamiltonians are 
[the variables /I(t) being the Dirac multipliers in front of the primary constraints 7Tx{t) ~ 0] 



'11 



i=\ 1=1 ^"^i i=l 

HD = T.77^-^-'^+ + f^-^>^- (2-3) 

i = l ^"^2 

The time constancy of the primary constraints implies the following secondary constraints 

fcxit) = {nxit),Hn} = K+^0, (2.4) 
which is the non-relativistic rest frame condition. 
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There are two first class constraints tF^ ~ 0, k+ ^ : A(t) and a center-of-mass variable 
are gauge variables. The Hamilton and Euler- Lagrange equations are [= means evaluated 
on the trajectories which minimize the action principle] 



rrii 



n^{t) = K+ ^ 0, 



N 



(2.5) 



2 = 1 



This is the non-relativistic limit of the relativistic rest-frame instant form of dynamics: 
Minkowski spacetime is replaced by Galilei spacetime and the Wigner hyperplanes are re- 
placed by the inertial observers seeing the isolated system is istantaneously at rest in the 
t = const, hyperplanes. 

Defining the non-relativistic center of mass 



l^ — Vi, 



(2.6) 



with m = Y^iLi f^i-i the gauge fixing qnr ~ implies A(r) and the decoupling of the 
center of mass, see Eq.(|2.9|). Instead, the gauge fixing r/+ = Y^f=i ~ does not imply 

A(r) ^ and the decoupling, just as it happens in the relativistic case f^. 

In analogy with the relativistic case of Ref. [^, let us introduce the following family of 
non-relativistic point canonical transformations 





pa 







defined by 



N-l 



—* 

Qnr 


Pqa = Pa 




T^qa 



(2.7) 



N-l 



V+ + -7ff 7aiPa = qnr + -J= ^ V aipqa, 

V A/ a=l ViV a=l 



They can be used inside the Lagrangian Ld\ the first one is the non-relativistic analogue of that 
used in Ref. Q. 

^^The total angular momentum of the N-body system is J = X^i^i 'Ti x = Qnr x k+ + Sg = 
ff^ X + S; S is the barycentric angular momentum or spin. 
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Ki = —K+ + VN > 'JaiT^a = + V iV > 7aj71",a, 

a=l a=l 



1 « 



TV 

j=i 

AT 

pa = VN^JaiVi, 



If ^ ^ ^ 1 /f "» ^- 



7r„ 



i=l 



Pqa Pa) 
1 

N-l N-1 N-l N 
S = Ypa^T^a= Y Pa^ T^qa + —?= Yl lakjPa X /«H 



a=l a=l 
N-l 

Sq = Pa X TTga, 

a=l 



^ mfc 1 ^ 

t^l^ ^ fc = 1 

TV N 
j=l j=l 

1 „ . m, 



E 7ai7aj = 5ij - — , E ^aiFai = 5ij . (2.8) 

a=l a=l 

Here, the 7aj's [and the Taj's] are numerical parameters depending on ]^{N — 1){N — 2) 
free parameters ||TT| , |2l|| . From now on we shall use the notation pqa for pa- 

Then, by using the equations of motion m[q^^.{t) + A(t)] = 0, we get the Lagrangian 
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Lrei and the Hamiltonian Hrei describing the relative motions after the separation of the 
center-of-mass motion 



N 



777 . r , 

i=l ^ 



1 



a=l 



2 o 



1 



1..N-1 



N 



kab 1 Tcj] — ^ba i Tc 



i=l 



1 ^ 



Eh=i mhlah Efc=i mk'-fbk 



N 



milailbi - 



m 



N 



Kb[m.,T,4 = NY: 



'~iai'~ibi 



i=l 

N 



N N -p N 
-L ai 



N -p -p 

rrii 



N 



k=i i=i k=i 



k=l 1=1 



i=l 

rrii 



^ N ^ N N 

+ 1^(E — ) E Tafc E Tfefc, 
i=l "'j /i=l fc=l 



4 

(t) = ^ kab[mi,r^i]pqi^{t), 



6=1 



L.AT-l 



re/ 



ab 



(2.9) 



The same result can be obtained by adding the gauge fixings g^r ~ which imply X{t) = 0, 
and by going to Dirac brackets with respect to the second class constraints k+ ~ 0, g^r ~ 0. 
The (6N-6)-dimensional reduced phase space is now spanned by pqa, T^qa = 7?a and from 
Eq.(|3) we have S = Sq = Efji' Pa x vf,,. 

At the non-relativistic level [y,^] the next problem for each N is to diagonalize the matrix 
kab[mi,Tai]- The off-diagonal terms of the matrix kab[mi,Tai] are called mass polarization 
terms, while its eigenvalues are the reduced masses. 



B. Jacobi normal relative coordinates. 



There is a discrete set of point transformations 

Pa - 



E ClbVgft 



(2.10) 



b=l 



which defines the relative Jacobi normal coordinates |^ and diagonalizes the kinetic energy 
term of the Lagrangian 
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k = 1, .., k]\r = ^{N — 1){N — 2); for two values ki and /c2 one has s^^^ = J2bLi^ ^, 



^N-l ri(kik2)4k2) 
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kab — kabirriijTai] — Cac^'^ (f^i^''^cd)Cdb\ 



cd 



Sea l^abSh 



4 



WW = i E f^i'^'^a^\t) = J 5: 57 



N-l 



Sk)2 



a=l 



a=l 



(2.11) 



where fj.^^^ are the reduced masses of the clusters and the ^ are called mass-weighted 
Jacobi coordinates. This form of the Lagrangian defines an Euclidean metric on the relative 
configuration space. 

The general Jacobi coordinates or vectors s^^-* organize the particles into a "hierarchy of 
clusters", in which each cluster, of mass iJ,'^\ consists of one or more particles and where 
each Jacobi vector s^^^ joins the centers of mass of two clusters, thereby creating a larger 
cluster; the discrete set {k = 1, .., k^) of choices of Jacobi vectors corresponds to the possible 
different clusterings of N particles. Usually, by "standard Jacobi coordinates" Sa one means 
the special set 



Sl 



S2 



\m1m2 



Vi - V2 



Ai(12) 



m 



/^(12) 

Nm 



N-l 



qai 



a=l 



mM2M3 
At((i2)3) \mif]i + m2f]2 



m 



/^((12)3) 

Nm 



L M2 

-rriiTai + m2r, 



a2 



a=l 



Mo 



-r 



aS 



SN-1 



rriN 



N-l 



[ TTliffi - Mn-iVN 



tuMn-iMn L ^ 

/^(...(12)3)...)Af-l) \l^i=l ^iV 



m 



M. 



- Vn 



N-l 



jj{kik2)T _ £){kik2) 1^ with the set of matrices d('^i'^2) (democracy transformations or kinematical 
rotations) forming the democracy group., which is a subgroup of O(N-l) [9,22| 

^'^ Mi = mi + Mi-i, Mn = m = J2f=i "^i! ^or the reduced masses we have /X(i2) = , At{(i2)3) = 
M3 
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N-1 Y^Af-l^r 

/^(...(12)3)...)jV-l) [ l^^=l mi^ai pi-. fr,.r,^ 

i M^_, ~ ^H/^'^- ^2.12) 

in which si joins particles 1 and 2 and is directed towards 1, while Sa is directed from the 
(a-l-l)th particle to the center of mass of the first a particles 0. 

Let us remark that our family [there are |(A^ — 1)(A^ — 2) free parameters inside the 7aj's] 
of point canonical transformations of Eqs. ( |2.7| ) contains as a special case the transition to 
the normal Jacobi coordinates of Eqs. (|2.13|) . 

An induced set of canonical transformations from the canonical basis Pqa, i^qa to the 
Jacobi bases is the following 

N-l N N-1 

^a'^ = E CSV-,. = v^E E cSV..^i, 

b=l i=l b=l 

N~l 1 N N-1 



-^^^=E-.^C^r' = ^EEr.dr" 

b=i ViV i=i b=i 



N-1 N-1 

"Sg = E PiJi^ ^ ~ E ^ '^sa\ 

a=l a=l 

N-1 ^(fc)2/ ^ 

a=l //ia 



]^ 1..N-1 



(k)-l^k) 
ah ' 



ab 



m 1..N-1 

^^^ + v^ E 7..d?vfi?. (2.13) 

ab 



C. More about the static orientation-shape SO(3) principal bundle approach. 

As said in the Introduction, the attempt of decoupling absolute global rotations from 
vibrational degrees of freedom led to the development of the static theory of the orientation- 
shape SO (3) principal bundle 0, which generalizes the traditional concept of body frame 
of rigid bodies. This bundle is a non-trivial (for > 4) principal S0(3)-bundle with the 
(3N-6)-dimensional shape manifold (with coordinates q^) as base and standard fiber SO (3) 
(parametrized e.g. by the Euler angles 6'°). For each given shape we need: i) the assignement 



^^See Ref. for more details on the general Jacobi coordinates and on special classes of them 
(like the Radau ones) treating the particles either in a more symmetric way or according to a more 
complex patterns of clustering; they are connected to the standard ones by kinematic rotations 
belonging to the democracy group. 
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of an arbitrary it reference frame; ii) the assignement of a body frame, identified by the value 
of the orientational variables 6'^ with respect to the reference frame. Recall that the 6'"'s 
are gauge variables in this approach. . 

A convention about which is the body frame for generic configurations of the N-body 
system, namely a local cross section of the non-trivial orient at ion- shape bundle, is equivalent 
to two independent statements: i) a choice of a set of S0(3)-scalar shape variables q^, 
yU = 1, .., 3iV — 6; ii) a choice of the explicit form of the components PaiQ^) [s^a^'^{q^''^^)] of the 
relative coordinate vectors with respect to the chosen body frame axes These components 
are connected to the coordinates pj^ in the space frame by 

The same relation holds for the components of every vector, like = R^'^{6°')Sg. If fr are 
the space frame axes and Cr the axes of the chosen body frame, we have Pga = P^afr = P^qa^r 
and Sq = Sgfr = SgCr for each possible shape of the N-body system. In particular the S0(3)- 
scalars have the same functional form in both space and body frames: A^B^ = A^B^ . See 
Refs. 0,^^-^ for the mathematical and physical aspects of the orientation- shape principal 
bundle approach. 

The main result from the theory of the orientation-shape bundle is that the transitions 
among different body frame conventions are interpreted as gauge transformations among the 
local cross sections of the principal bundle. Therefore, a gauge transformation is a shape- 
dependent proper rotation S{q) G 5*0(3) that maps the body frame with axes ir into the 
body frame with axes 0. Instead of this passive change of coordinates on the fibers, 
one can consider an active (gauge dependent) right action of S0(3): {R,q^) ^ {RQ,q^), 
Q e 5*0(3). The corresponding symplectic right action in phase space, associated with the 
left-invariant vector fields on SO (3), is generated by the non-conserved body frame spin 
components 5"^. On the other hand, the left action of S0(3) p^^ ^— Q^^ Pqai {^il^) ^ 
{QR,q'^), Q G 5*0(3) is generated by the space frame spin components 5^ (Noether 
constants of motion), associated to the right-invariant vector fields on SO (3). As already 
said, it holds {5g^ S'q} = e-"5^^ {S^, S^} = -e"-5^", {5^, S^^} = 0. 

In conclusion, within the static orientation-shape bundle approach the 6A^ — 6 shape 
variables q'^ are gauge invariant quantities, because their definition does not depend on 
the body frame convention, while the body frame components of any vector are gauge 



^This is the choice of a gauge for the orientation variables, independent of the shape coordinates 
q^; for each shape, one gives the positions of the N particles relative to the body frame axes e^. 
The orientation variables, for example the Euler angles 0", identify an S0(3)-element R{9") in the 
fiber over the given shape q^; the reference orientation for each shape is such that /)^^ = /jg^. 

20We have = Srs{q)e„ R = R'S^{q), A = A'-fr = i^'e^ = i'^'e^, A"- = R'-'A' = R'^-'A''. The 
new orientation angles 9 " depend on the old ones 6°^ and on the shape variables q'^ too. 

^^I.e. the action of the structure group on the bundle, which is independent of the choice of any 
cross section and is called a 'gauge-invariant action'. 
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quantities. On the other hand, let us stress that only the vectors independent of the body 
frame convention have their space frame components as observable physical quantities. The 
angular velocity a; is a clear instantiation of a body frame dependent vector, for which both 
the body frame cD^ and the space frame uf components are gauge quantities. Note that they 
are not even gauge covariantf^, because under a change of body frame R{6°') ^ R {O"") = 
it holds tu'' = S'''[iu'' - 7^g^] with e'^'^'-f'^q'' = [S^Sy, so that lu'' ^ J''. 
These results suggest to consider local point canonical transformations of the form [/i = 



l,..,3Ar-6] 



qa 



"a 



qa 







Pa 





(2.15) 



They define a canonical basis, in which the local "orientation" coordinates 6°', a = 1, 2, 3 
are either the Euler angles or any other parametrization of the group manifold of SO (3). 

D. Non-Relativistic Rotational Kinematics in Jacobi Coordinates 

In this Subsection, using Jacobi coordinates, we shall elucidate the Lagrangian and 
Hamiltonian treatment of the orientation-shape bundle approach. In Appendix B a reformu- 
lation of these results is given in terms of arbitrary (non- Jacobi normal) relative coordinates. 

Given a set of Jacobi coordinates s^'^^'', let us introduce the associated body frame coor- 
dinates s'^^^'^ and velocities v^''^'^ 



The body frame components of the angular velocity u are 

a;'-(r,r) = --e"™[i?^i?]™(^",r) = ie™^[i?^i2]™(r,^"). 



(2.16) 



(2.17) 



As said before, also the space frame angular velocity components uj'^{6°',6'^ 
Jl^s^Qa^ ^s^Qa^Qa^ are gauge dependent. 

The Jacobi momenta are 
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A gauge covariant quantity transforms as A^' = S'^'^{q)A 
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^{k)r _ (k) 



sa 5 



For the spin we have 



N-l 



a=l 



5,= E4' 



iV-1 



(fc) jC^) X 



a=l 
Af-1 

E/^i' 

a=l 



a=l 



(fc) 



By introducing the Euchdean tensor^ 

fuv (Jk]\ _ ^k)2s.uv _ ~{k)u^{k)v 
-'{aa)W J ~ ^ i 

and the body frame barycentric inertia tensor 



N-l 



a=l 



we get the following expression of the body frame spin components^ 



N-l 



a=l 
N-l 



a=l 

j{k)uv^^(k)^ ^{k)^ J^. ^ ^W) gW/^ 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



^^In the relativistic case where the tensor 
be defined. 

^■^This tensor is gauge covariant, while J^'^)" 
tensor I™ = j^^-i' j^"^ j'^s gauge invariant. 

25pQj. = we get the rigid body result 5^' 
no longer hold in the relativistic case 



f inertia does not exist, only the tensors ( p.20| ) can 
= S'^'''{q)i^''^''''''{S'^Y"{q); the space frame inertia 
= i^'^y^uj^ . Let us remark that also these relations 
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where 



N-l 

a=l a=l 

/iW) i(f)"(g('=), /i^'^)). (2.23) 



The quantity A^^\q^''\ fx^''^) is the S0(3) gauge potential of the orientation-shape bundle 
formulation 0. Note that it is not gauge covariant: A^f^'' = S'^^{q)[A^^'^''^ + 7^]. Its field 
strength (curvature form), called the Coriolis tensor, is the gauge covariant quantity [B^J'^ = 



m) = ^^^^^ _ _ m X m 



(2.24) 



Let us stress that Eq. (|2.22|) does not provide an effective separation of rotational and 
internal (vibrational) contributions to the angular momentum, since the separation is gauge 
dependent within this approach . 

Eq.( |2.22| ) can be inverted to express the body frame angular velocity components in 
terms of the body frame spin components and of the gauge potential 



(2.25) 



The non-relativistic Lagrangian for relative motions can then be rewritten in the following 
forms 



1 



l..Af-l 



n,)'Jk)2 



-'rel 



E, u u _ 1 ^v-^"^ (A;) Uk) _ ^ sr^ (k) ^vt', 

ah a=\ a=l 

1 



j{k)uv^^u ^ A{k)u^(k)^,■^|^^v ^ A{k)v^(k)u-^ ^ gllk^'^^^q^''^'' 



(2.26) 



where 



a=l 



^l^q{k)i^q(k)u^ 



(2.27) 



^^See Ref. for the monopole-like singularities of the gauge potential at the N-body collision 
configuration. 
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is a pseudo-metric on shape space 0, while 



9tKq^'\ f^^'^) = (2.28) 

is a true gauge invariant metric on shape space [g^^] = gjj'J'] 0. 

A manifestly gauge invariant separation between rotational and vibrational kinetic en- 
ergies is exhibited only in the last two lines of Eq.( p.26| ) 0. In order to clarify this point, 
a velocity multivector {6"',q'^) was introduced in Ref. together with its anholonomic ver- 
sion (u)'", g^) in which Euler angle velocities are replaced by the body frame angular velocity 
components. A metric tensor for these multivectors is naturally induced by the Euclidean 
metric of the kinetic energy. The intrinsic notion of vertical vector fields of the S0(3) prin- 
cipal orientation-shape bundle corresponds to the purely rotational velocity multivectors 
{u^, 0) defined by the gauge invariant condition q'^ = 0. On the other hand, there is no 
gauge invariant definition of purely vibrational velocity multivectors^, since any such defi- 
nition is connected to a horizontal cross section of the principal bundle and, therefore, to 
the assignement of a connection form. Each connection gives a definition of horizontality 
and the possibility, through a gauge fixing, to choose a certain horizontal cross section as a 
connection- dependent definition of vibration. 

In Ref. it is shown that a special connection C can be defined by requiring that the C- 
horizontal vector fields are orthogonal (in the sense of the multivector metric) to the vertical 
ones and that the associated C-horizontal cross sections (defined only locally for > 4) are 
identified by the vanishing of the body frame spin components S"^ = ^ By privileging the 
connection C, we get the following splitting of an arbitrary velocity multivector into vertical 
and C-horizontal parts 

(d;^ga = (^'^ + 4'^V,0), 
(u;^g^)c. = (-i?V,g^). (2.29) 



^'^It is neither gauge invariant nor gauge covariant: h^/S = h^^} +^^^1^^^'^^ a!i^ ^ +^'^1'^^^''^'^ A^^^ ^ + 

2^It can be shown [§ that the inverse metric is g^"' = Y.a=i " l^- 
29For = we get the rigid body result L = \i^''>'Cj''-uj'' = \I^^>'uj''oj'' . 
^'^For instance the simplest choice w'' = is clearly not gauge invariant. 

^^A system velocity is C-horizontal if and only if the associated spin vanishes and horizontal 
vector fields describe purely vibrational effects in a gauge invariant way; this also implies that 
the C-horizontal cross sections cannot be interpreted as (6A^ — 6)-dimensional submanifolds of the 
configuration space. 
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It is just this C-splitting which identifies the metric ( |2.28|) on shape space and the 
manifestly gauge invariant separation of the kinetic energy given in the last lines of Eq.j^^Wi). 

The standard fiber of the orientation-shape bundle is SO (3). Its group manifold admits 
many parametrizations. If one uses the local parametrization given by the Euler angles 6" 
Q the form of the right-invariant vector fields and 1-forms on SO (3) Q is 



d 



'OOP'' 



(a) 



OR 



<5l 



(2.30) 



If the Euler angles are defined by the convention R{9'^) = Rz{6^)Ry{6'^)Rz{0^), one has 



/ —sinO'^cosO'^ sin 6^ 

sin9^sin9^ cos 9^ | and det A^^^ = —sin9'^ 
2 1 



\ cos 9^ 



The left-invariant vector fields on SO (3) are {[xi^} , X^^^ 



(a) '^-(/3)J - 0) 



(2.31) 



The linear relation between the body frame angular velocity and the velocities 9" is 

cj^ = A(^^f='-\9^)9'^. (2.32) 
Using the anholonomic components {uj^, q^) Q of the velocities instead of the holonomic 



^^See Ref. ||2^ for the parametrization of the S0(3) group manifold with a 3-vector e determining 
the rotation axis and the rotation angle V' by \e\ = 2sin ^ . . 

^^Recall that they are the generators of the infinitesimal left translations on S0(3) and are thought 
as body frame quantities. 

^"^Recall that they are the generators of the infinitesimal right translations on SO (3) and inter- 
preted as space frame quantities. 

^^Using the dreibein given by the right-invariant vector fields on S0(3) and regarding the La- 
grangian as function of 6", cD'', q^, [see Eq.( |2.26| )]. 
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basis {6", q^), one gets the following canonical anholonomic momenta fj\ 

dLrel 



or 



(k) _ 9Lrel _ ik)-(k)u , r? . /(fe) 



^9^" 95,, dSgrde^^ ' ^dSg dSg 



'df_dg_ _ df dg 
dpf, dpf, dqi 



Ch, 



(§1, p,). = {Sl, S,-Af{q^'\^,^'^)), 

K P,)cH = {0, - 5^, -Ajng^^), (2.33) 



The last lines show the decomposition of the momenta into vertical and C-horizontal 
parts. 

Finally, the Hamiltonian becomes 

1 A^-l -(fc)2 

J^rel — r) 2^ (f,) — ^ P/i Q ^rel — 

^ a=\ /ia 

= ^[^g(/^'^^"'(g^^\yu(^^))™^; + 

+ [g^^'^^^vf - 5^, • Af){j>^:;^ - S, ■ 4^))) iq^'\f^^'^)] . (2.34) 

For q^ = 0, namely p^ = Sq- A^l^\ one gets the rigid body Hamiltonian for pure rotations 
without vibrations H = l{I-^Y''>''S'^S'g = |(/-i)('=)^'^S^5g^ 

See Appendix A for the form of the Lagrangian and Hamiltonian equations of motion in 
the orientation-shape bundle approach and Appendix B for the reformulation of the results 
of this Subsection with arbitrary coordinates. 



^^The body frame spin components replace the momenta vr^Q = A^^^i^ Sqr conjugate to 9". 

^^Let us remark that one could also use anholonomic gauge invariant shape momenta = Pi_i — 
Sq ■ with Poisson brackets: = {Sq,p^} = x Sq, = <5^^, 

{p^l,Pu} = Sq • Bj^J. 
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III. CANONICAL SPIN BASES. 



The static theory of the orientation-shape SO (3) principal bundle is based on point 
canonical transformations of the type of Eg. ( p. 151) . 
Following the preliminary work of Ref. 



]15|] , we look for a set of non-point canonical 



transformations from the relative canonical variables pqa = pa = \ N X]i=i la 
-^Y.i^=i^aii^i of Eq. (|2.8| ) to a canonical basis adapted to the S0(3) subgroup of the 
extended Galilei group and containing one of its invariants, namely the modulus of the spin. 
Jacobi coordinates will not be used in this Section and the comparison with the orientation- 
shape formalism has to be done by using Appendix B. 

Again, the configurations with 5* 7^ and with 5 = have to be treated separately. The 
special connection C of the static orientation-shape SO (3) principal bundle is not included 
in our description, which is valid only for the S configurations. Accordingly, after the 
exceptional case N=2, we shall find that in the case N=3 the results of the static trivial 
orientation-shape SO (3) principal bundle are recovered in a gauge of the xxz type. On the 
other hand, our results for > 4 will differ substantially from the static non-trivial SO (3) 
principal bundle approach. 



A. 2-Body Systems. 



The relative variables are pq = p, ifq and the Hamiltonian is Hrei 

is the reduced mass. The spin is Sq = pq x vf^ [Sg = 
Let us define the following decomposition PI 



Pq = PqR, 



Pq 



Pq 



Pq, 



TTqR - 



Sq 
Pq 



Rx 



_ 5" - 

^qPq -pq ^ ^q^ 

Pq 



TVq ■ R 



T^q ■ Pq, 



Q — Zl 



where p 



mim2 
mi+m2 



R' 



S„-R = 0. 



(3.1) 



Therefore, besides the standard space or laboratory frame with unit vectors fr, we can 
build a spin frame, whose basis unit vectors Sg, R, Sg x R are identified by the 2-body 
system itself. Since {S^, S^} = e'^'^S^, {R\W} = 0, {R\ S^} = e^^'^R^, Sq and R are the 
generators of an E(3) group containing SO (3) as a subgroup. The E(3) invariants turn out 
to have the fixed values R^ = 1 and Sq ■ R = 0. 

Usually one considers the following local point canonical transformation [like in 
Eqs. (|2.15| )j to polar coordinates 



pq 



9 if 


pq 


T^qd T^qf 


T^qp = TTq 



(3.2) 



^^The notation R for the unit vector pg is used for comparison with Ref. |15| 
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Pi 
Pi 
p\ 



PqSin 6 cos if, 
PqSin Osin 

PqCOS 9, 



TC 



r.1 



sin OCOS ^T^qp — 



sm Lp 
PqSin 9 



COS ipcos 9 



-7T„ 



COS (f sin ipcos 9 

Sin9sinip'Kqp + ^ ITq^ H 7r„ 



PqSin9 



cos 9tt 



sin 9 



qp 



-IT. 



7r„ 



2 1 



qe, 



TC 



qip 



sin^ 9 



cos 9cos if 

sin 9 
cos 9sin if 



TXq^ - SmLfTlgg, 



Sin I 



St = TT, 



9¥" 



Si 



TTq^ + COS ipHqQ, 



2 , ^gy 
'^'^ stn^9' 



Pq 



C2 



(3.3) 



pq 



'P^q, 



n Pq 

COS 9 = — , 

Pq 

pI 

t9^ = —, 

pq 



sin 9 



simp 



pq 
Pi 

^{Pq? - {pW 



tg9 



cos if 



^{Pq? - [Pl? 

Pi ' 

^{Pq? - {pD 



"^qp '^q ) 



TTqd = Pq[cos 9{cos ipTT^ + sinifTi'l) — sm 6*7?^] 



^|{Pqf^W' 



TT, 



PqSin9{cos^7il - sin^n]) = pln^ - pjvrj = Sl 



(3.4) 



After this point canonical transformation, the Hamiltonian becomes Hrei = ^ 
while the static shape canonical variables are the pair pq, iTqp. 

As shown in Ref. instead of this point canonical transformation, it is instrumental 
to consider the following non-point canonical transformation adapted to the SO (3) group 
valid when Sq ^ 



^^Note that in the new canonical basis the invariant Sq = \Sq\ becomes one of the new canonical 
variables. 
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7r„ 



a (3 


Pq 


Sq 


^1 



(3.5) 



where 



a 



-1 1 
t9-^[pq- 



7r„ 



(P.)^3 



3 



1 s „ 

/3 = tr'^, 



C2 



COS /3 



91 



(3.6) 



The two pairs of canonical variables a, Sq, (3, S^ form the irreducible kernel of the 
scheme ^4 of a (non- irreducible, type 3, see Ref. ITH]) canonical realization of the group E(3) 



generated by 5*^, R, with fixed values of the invariants R^ = 1, R ■ Sq = 0, just as the 
variables S^, (3 and Sg form the scheme A of the S0(3) group with invariant Sq. 
Geometrically we have: 

i) the angle a is the angle between the plane determined by 5*^ and /a and the plane 
determined by 5*^ and R] 

ii) the angle (3 is the angle between the plane Sq - /s and the plane /s - /i. 



We have 














= pyq 


- py = 


— sirnpTTqg — COS ipCOtg OTTg^ = 


^{Sqf 


-{SlfcosP, 


Sg 


= pyq 
= py 


- py = 

- py = 


COS ifKqe — sin Lpcotg Oj^q^ = y 

T^qip-i 


'{Sqf - 


{SlYsmP, 


Sq 




sin? 6 









(3.7) 



R^ 


= P\ = 


R' 


= p\ = 


R' 


= P\ = 


ky 


OgtX 


Rf 


= Sglx 



{Sq X RY 



Sq 

sj 



c-3 

s: 



q V \-qy y-qy 
Oq 



S^„ 



{Sq X R) = SgR — SgR = sin f3cosa + -^cos Psina, 
{Sq xRY = S'-'R' - SlR^ = -cos (3cos a + -f szn (3 sin a, 



Sa 



-SIR' - S'gR' = -^^{SqY - {Slfsma, 



- , dR{a) - , TT, 

Sq X R{a) = = R{a + -), 
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^ a = -tg-' . ^ . (3.8) 

From the last line of this equation we see that the angle a can be expressed in terms of 
and R. Given the Hamiltonian description of any isolated system (a deformable body) in its 
rest frame with conserved spin S{q, p) [q, p denote a canonical basis for the system], a solution 
a{q,p) of the equations {a{q,p), S{q,p)} = 1, {a{q,p), (3{q,p)} = {a{q,p), S^{q,p)} = 0, 
allows to construct the unit vector R associated with the isolated system and then to build 
the spin frame and the E(3) group. 

The following inverse canonical transformation holds true 



Pq = PqR{(^,P,Sg,Sgj, 

, - n.Ria, (3, S,, 5j) + S,, S'^) x R{a, (3, S,, 



TT, 



Pq 



2 



^ ^q = K + §• (3-9) 

Pq 

In this degenerate case, the dynamical shape variables pq, 7fg coincide with the static ones 
and describe the vibration of the dipole. 

The rest-frame Hamiltonian for the relative motion becomes [/ is the barycentric tensor 
of inertia of the dipole] 



Hrcl — ~ 

2 



lr^_i^o TT^T , 2 mim2 



PP . P= : , (3.10 

nil + ^2 



while the body frame angular velocity is 



dHrel Sq /n tt\ 

uj = — z — = 3.11 

dSq I ^ ' 

We conclude this Subsection with some more details on what has already been anticipated 
in the Introduction concerning the canonical reduction. When Sq ^ 0, Eq.( |3.5| ) explicitly 
shows that a non-Abelian symmetry group like SO (3) does not allow a canonical reduction 
like in the Abelian case of translations. In this latter case we can eliminate the three Abelian 
constants of motion k+ ^ and gauge fix the three conjugate variables qnr ~ 0. In the 
non-Abehan case we could surely fix S^, [3 by imposing second class constraints 5^ — a ~ 0, 
/3 — 6 ^ 0, and eliminate this pair of canonical variables by going to Dirac brackets. However, 
since a is not a constant of motion [see later on Eq. (|3.1CI|) ; we get instead ^{Sq — 5") =0], 
we can only add by hand the first class constraint 5*^ — 5 ~ (5* 7^ 0). It is only after the 
solution of the equations of motion, that we could also complete the reduction by adding 
a - oisoiution ^ as a gauge-fixing. 

In absence of interactions the solution for a can be easily worked out. The Hamil- 
ton equations, equivalent to Pg = 0, are f>q = 'nq/n, vfg = (S'q)^//ipg, 5*^ = 0, a = Sq/ ppg. The 

solution pq = bt + a [a, b constant vectors] implies: pq = \bt + a\, R= {bt + a)/\bt + a\, 
^q = 2pb- {bt + a)/\bt + a\, Sq = pV2\axb\,a = arccos [k^Sq/^iSq^ - (Siy). 
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In the 2-body case the condition Sg ~ 0, imposed as three first class constraints, is 
equivalent to Pg — kifg ^ and selects only the solution pg{t) = Ae^^ [k, A and B are 
constants]. 

B. 3-Body Systems. 

In the case N=3 the range of the indices is z = 1,2,3, a = 1,2. The spin is Sg = 
Z]a=i Pqa X T^qa = J2a=i Sga after the canonical transformation which separates the internal 
center of mass 





pqa 




T^qa 



(3.12) 



The relative motions are governed by the Hamiltonian 

1,2 



1 



cil'^qa ' T^qb- 



(3.13) 



a,b 



Again, we shall assume Sg ^ 0, because the exceptional SO (3) orbit Sg = Q has to be 
studied separately. This is done by adding S'g ~ as a first class constraint and studying 
the following two disjoint strata with a different number of first class constraints separately: 
a) Sg ~ 0, but Sgi = —Sg2 7^ 0; b) Sga ~ 0, a = 1, 2 [in this case we have pga — kaT^ga - 



For each value of a = 1, 2, we consider the non-point canonical transformation (|3.5| ) 



Pqa 



qa 



Oia (^a 


pqa 


C C3 
"^qa ^ga 


TTqa 



(3.14) 



where 



Pa = tg 



s 



qa 



C2 
-1 '-'qa 

CI ' 
>-'qa 



Pqa 



(p. 



qa 



sin Pa 



qa J 



Pqa 



-vr 



qa 



q2 



{s,ar-{siar 



cos Pa 



{s,ay 



(Sir 



(3.15) 



Let us remark that with the canonical basis (3.14) the degenerate case defined by imposing the 
constraints Sg ^ with Sgi ^ —Sg2 ^ implies the three extra first class constraints Sgi — Sg2 ~ 0, 
Sgi + Sg2 ~ 0, /3i — /32 ~ 0: therefore we get three arbitrary conjugate gauge variables ^(ai — 02), 
^(A + /32), ^5*^1 + 3^2) / 0- Besides these three pairs of conjugate variables, a canonical basis 
adapted to « (with Sgi —5*^2 / 0) also contains the physical variables a = 01 + 02, 
S = i(5gi + Sg2), Pga, T^qa, and the Hamiltonian (|3.13| ) becomes Hd = ^Y^alKb^^qa ' ^'?fels,=o + 
Ai(t)(5gi - Sg2) + X2{t)iS^i + 5^2) + A3(t)(/3i - P2) [the A's are Dirac multipliers]. 
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p / ^ p Pqa - p2 1 

P(ja Pqa-l^ay Pqa y Pqai ^a Pqa) J^a 

Pqa 

T^qa = T^qaRa + —Sqa y< Ra, T^qa = T^qa ' Ra- (3.16) 

pqa 



Pqa PqaPqa{o!a) Pai Sqa, S^^^ PqaRa{cta, Pai Sqa, S^^^, 

3 "S* a " 3 3 

T^qa T^qaPqaic^a^ Pay Sqa, Sg^) H Sqa{Pa, Sqa, S^^ X Pqaiflia, Pa, Sqa, S^^)) 

Pqa 

~ " 1, S a " 3" 3 

= T^qaRaiS^a, Pa, Sqa, S^^) H Sqa{Pa, Sqa, S^^) X Ra{oia, Pa, Sqa, S^^). (3.17) 

Pqa 

We have now two unit vectors Ra and two E(3) realizations generated by Sqa, Ra respec- 
tively and fixed invariants R^ = 1, Sqa ■ Ra = (non-irreducible, type 2, see Ref. |T5||). 

Then, the simplest choice, within the existing arbitrariness (footnote 10), for the or- 
thonormal vectors and x functions only of the relative coordinates is 



N = liRl + R2) = liPql + Pq2), N=^ \N\ = ^^^Pf^, 



X = l{Rl- R2) = liPq^ - Pq2), X=^y IXI = y ^-^f^ = 



^ X X = -2P91 X Pg2, |iV X xl = |iV||x| = -/l^"(^^7^7^, 

iv ■ X = 0, {iv^ N^} = {x^ x'} = {iv^ x^} = o, 

Rl=Pql=N + X, R2=Pq2=N-X, Rl ' R2 = Pql ' Pq2 = - f ■ (3.18) 

Likewise, we have for the spins 



Sq 


— Sql + 5*^2, 




Wq 


= Sql — Sq2, 




— * 

Sql 


= liSq + Wq), 


— * 

Sq2 






rsu QU 



1 

5' 



(3,19) 



We therefore succeeded in constructing an orthonormal triad (the dynamical body frame) 
and two E(3) realizations (non-irreducible, type 3, see Ref. Il5|): one with generators Sg, N 



41c 



See Ref. |15| with the interchange pga ^ TTqa in the canonical transformation introduced there. 
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— # — * ^ — * 

and non-fixed invariants |A^| and S ■ N, the other with generators Sg and x ^i-nd non-fixed 
invariants \x\ and 5*^ ■ x- As said in the Introduction, Eg . ( |1 . 1|) , this is equivalent to the 
determination of the non-conserved generators of a Hamiltonian right action of SO (3): 

Si = Sg ■ X = Sg ■ = Sg ■ N X X = Sg ■ 62, = Sg ■ N = Sg ■ 63. 

The realization of the E(3) group with generators Sg, N and non-fixed invariants A^^, 
Sg ■ N leads to the final canonical transformation introduced in Ref. [|T5| 



where 



Pqa 



qa 



"1 Pi 


"2 P2 


pqa 


Sql 


Sg2 S^2 


T^qa 





\N\ 


Pqa 


Sg = Sg Sg Sg = Sg ■ N 




T^qa 



(3.20) 



'1 ■ Pq2 



Sg — Sg ■ N — /- Pga 

a=l 



X TT, 



COS i/j = Sg ■ N 



SI 

S-' 



qa 



simp 



Pql + Pq2 
1 + Pgl • Pq2 
1 

s 



SgCOS 



{Sgf 



Sn S, 



s' 



Y.pi'^ 



X 71, 



qa\, 



a=l 



II(Pga X IT. 



3 

qaj , 



a=l 



a 



-tg 
-tg 



-1 



{Sg X Nf 



[Sg X {Sg X N)f 
-1 [Sq X {Pgl + Pg2)Y 



[Sg X {Sg X [p^l+^^2])]3' 



s^ 



7 = tg 



-1 Sg-{Nx x) 



Sq ■ X 



tg 



02 

1 

CI ' 

Dg 



smy 



tg-' 



C2 



(C3^2- 



{SgY 



COS 7 



g ■ Pg2 X Pql 



1 + Pgl ■ Pg2 Sg ■ {Pgl - Pg2) 

Wg-{Nxx) _Wg-{Nx x) 



\x\ 



1 - 



V2ELl(-)"+Vga X ^ga ■ {Pg2 X Pgi) 



[I- pql - Pq2]J^+ pql- Pq2 



ql 
^q 



{Sq? 



{s'y 



(3.21) 
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For N=3 the dynamical shape variables, functions of the relative coordinates Pqa only, 
are |A^| and pqa, while the conjugate shape momenta are ^, Tiqa. 

The final array (|3.2CI| ) is nothing else than a scheme 5 of a realization of the E(3) 



group with generators Sq, N (non-irreducible type 3). In particular, the two canonical pairs 
Sg, /?, Sq, a, constitute the irreducible kernel of the E(3) scheme A, whose invariants are S^, 

\N\; 7 and ^ are the so-called supplementary variables conjugated to the invariants; finally, 
the two pairs Pqa, liqa are so-called inessential variables. Let us remark that S^, P, Sq, a, 7, 
^, are a local coordinatization of every E(3) coadjoint orbit with 5*^ = const., \N\ = const. 
and fixed values of the inessential variables, present in the 3-body phase space. 

We can now reconstruct Sq and define a new unit vector R orthogonal to Sq by adopting 
the prescription of Eq. ( p.8|) as follows 



sl = ^^{Sqy-isiycosp, 

Oq 



S'q = ^^{Sqy-{Siysznf3, 



Sq 
Sq 



1 s ^ 

R =sin(3sina — ^cos(3cosa, 

Oq 

R = —cos j3sin a — -^sin (3cos a, 

Oq 

1 



Sq 



R' = -x/iSq)^-iSiycosa, 



R' = l, R-Sq = 0, {R'',R'} = 0, 



5-3 

{Sq X R)^ = SgR^ - S^R"^ = sin (3cos a + -^cos (3 sin a. 



S^ 

{Sq X i?)^ = SgR^ — SgR^ = —cos (3cos a + -^sin (3sin a, 

Oq 

1 



{Sq X Rf = - = —^{SqY - {S^ysina, (3.22) 

Oq 

The vectors Sq, R, Sq x R build up the spin frame for N=3. The angle a conjugate to 
5'g is explicitly given byQ 

1 (So X iV)3 . {Sq xRf , , 

[SqX{SqXN)f [SqX{SqXR)f 



The two expressions of a given here are consistent with the fact that Sq, R and N are coplanar, 
so that R and differ only by a term in Sq. 
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^45 a consequence of this definition of R, we get the following expressions for the dynam- 
ical body frame TV, x, TV x % in terms of the final canonical variables 

TV = cosipSg + sini^R - ^.3^ + ^^^3^)^ - ^s^^R = 



X — sin ipcos ^ySq — cos ipcos jR + sin jSq x R — 

I - 

= Trv('S'q)^- {Sqycos'ySq- -^cos'yR + sin'ySq x R 

^q ^q 
_ ^ 3 _ Sq R + S'^ Sq X R ^ 
~ Sq Sq 1 (^3)2 " 

= x[Sq,a;Sg,l3;Sg,-f], 



N X X — sin ipsin 'jSg — cos ipsin jR — cos 'jSq x R — 

= ^\/{Sqy-{Siys^rijSq - ^sinjR - cos^Sq X R 



Sq Sg Sq Sq X 



- 

= {Nxx)[Sq,a;Sl,P;Sl^], 

— sin i/jcos 7X + sin ipsin '-/N x % + cos ipN 



def 1 



Six + S'N X X + S'N 



R — —cos i/jcos 7X — cos ipsin 7TV x x + sin ipN, 



Rx Sq= -sin 7X + cos 7TV x x- (3.24) 

While ip is the angle between 5"^ and TV, 7 is the angle between the plane TV — x and the 
plane Sq — TV. As in the case N=2, a is the angle between the plane Sq — /s and the plane 
Sq — R, while (3 is the angle between the plane 5'^ — f^ and the plane fa — fi- See the Figure. 
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FIGURES 




Frames 

FIG. 1. Space frame (/r), spin frame (S*, R, S x R) and dynamical body frame {N, x, N x x) for N > 3 
particles. V - angle between N and S. a - angle between the planes S-f^ and S-R. fi - angle between /i 
and the projection of S onto the plane /1-/2. 7 - angle between x and the projection of S onto the plane 



Owing to the results of Appendix C, which allow to re-express Sqa = \Sga\, S^^, Pa = 

52 

tg~ ^ in terms of the final variables and owing to Eqs. p.lq ), ( p.23|) which allow to get 



"a 



we can reconstruct the inverse canonical transformation. 



{Sqa X {Sqa X Ra ) )^ 

The existence of the spin frame and of the dynamical body frame allows to define two 
decompositions of the relative variables, which make explicit the inverse canonical transfor- 
mation. For the relative coordinates we get from Eqs. ( |3.18| ) and 



= [Pga ■ Sg]Sg + [Pga " R]R + [Pga " Sg X R] S g XR = 

_ pga 



Cl c3 



R- 



\a+l 



02 

^q 



-.S„ X R 



Pga[Sg,a;Sg,p;Sg,-f;pga, \N\]. 



(3.25) 
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The analogous formulae for the relative momenta are [see Eg. (|C5|) for the expression of 
the body frame components of vfga] 

T^qa — '^qa-tia H t^qa X rCa — T^qaPqa H Jqa ^ Pa — 

Pqa Pqa 

= [n,a ■ N]N + [nga ■ x]x + l^qa ■ N X x]N X X = 



Ka ■ S,]Sg + [nga ■ R]R + l^ga ■ Sg X R]Sg X R = 

^ [{{TTga ■ N)Sl + (7f,, ■ x)Sl + (vf,, ■ N X x)Sl)Sg + 



+ {{TTga ■ iV)v/(5,)2-(S3)2- [(vf,„ . x + {^ga ' N X x)] 
+ [{T^qa ■ X) - (7?9a ' N X x)] 



Cl C3 



(Sg)' 



R + 



C2 c3 



-.R X S„ 



iSq^' 



'qj V-g. 

= TTqa[Sq,a;S^,p;S^,-f; |iV|, ^; p^^, vf^J. (3.26) 

Finally, the results of Appendix D allow to perform a sequence of a canonical transfor- 
mation to Euler angles a, f3, 7 with their conjugate momenta, followed by a transition to 
the anholonomic basis used in the orientation-shape bundle approach H 



a 



(3 7 



Q C3 C3 
^q >->q >Jq 



a /5 7 



Pa Pp Pi 



non can. 



a /3 7 



cl C2 c3 
'-'q 
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7 



Sq = ^/iSlV + (S^g)' + (53)2, 

—sin Pcos 'jS^ + sin (3sin jS^ + cos PS^, 
Pptg(3 



arctg 



C _ P&Pj ' 

9 SgCOSP 



IT 

2 



7 — arctg 



ctg /3p, 



sin (3 



P = a + arctg 



ctgPPa-£^ 



71 

2' 



(3.27) 



Here pa, Pp, are the functions of a, /3, 7, S^ given in Eqs.(^). The equations (p3p, 
( p. 27] ), ( p.l8| ) and S"^ = Sg ■ N x x lead to the determination of the dynamical orientation 
variables a, P, 7 in terms of pga, Tiga- Let us stress that, while in the orientation-shape 
bundle approach the orientation variables 9" are gauge variables, the Euler angles a, /?, 7 
are uniquely determined in terms of the original configurations and momenta. 

In conclusion, the complete transition to the anholonomic basis used in the static theory 
of the orientation-shape bundle is 



a 7 


\N\ 


pqa 


C c c3 c3 
i->g — iJq iJq iJq 




T^qa 



non can. 



a /? 7 


\N\ 


pqa 


cl 62 63 
Dq Dq 




T^qa 



(3.28) 
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In order to further the comparison with the orientation-shape bundle approach, let us 
note the following relation between the space and body components of the relative coordi- 
nates. Eqs.(g), ( p:28|) , ( PH) and (^) imply 



P\a{(l) 



l-iV2, 



and 

so that the final visualization of our sequence of transformations is 



(3.29) 
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re, 



qa 



non can. 



a /3 7 


q^Pqa) 


cl 62 c3 
'-'q '-'(7 


Pt,{pqa,Trqa) 



(3.30) 



Note furthermore that we get p^^ = Pqa ■ x % = by construction and this entails that 
using our dynamical body frame is equivalent to a convention {xxzz gauge) about the body 
frame of the type of xxz and similar gauges quoted in Ref. 

Finally, we can give the expression of the Hamiltonian for relative motions^ in terms of 
the anholonomic Darboux basis of Eqs.( p.2'^ ). By using Eq. (|C7|) we get 
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PqlPq2 
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(5f^ 
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1 



2^ Pi 

,k-^ 



k'^ 



22 > 



Pg2 



PqlPq2 



+ 2k-\2\N\^l-N^{'-^ 



Pq2 



Pql' 



11 



2\N 



N2 Pql 



2 )^q^q + 
Pq2 



4p^i 



22 



'^q2 + 



4p^2 



+ 



These gauges utilize the 'static' shape variables pqa, 4> with cos (f) = 2N^ - 1, which can be 
expressed in terms of the dynamical shape variables pga, 1-^1- 

^^The Hamiltonian in the basis (|3.20| ) can be obtained with the following replacements Sg = 
y^(5,)2-(53)2cos7 and = ^ {S^y - {S^q)^sinj. 
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+ 2k 



-1 



12 



(2iV^ - l)vf,ivf,2 - |iV|(l - N')iC^ + ^) + 

Pql 



Pq2 



+ 



4PglPg2 



def 



def 1 

~ 2 



9 



fJ,U 




(3.31) 



where = {pgi, Pq2, p^ = {'rcqi,'rcq2,0 t^ie dynamical shape variables. In the last 
two lines we have rewritten the Hamiltonian in the form of Eq. ( 



In Appendix E we evaluate the quantities A^^i^q), g^'^iq), i~^^'^{q) appearing in the stan- 
dard static theory of the orientation-shape bundle in the xxzz gauge, adopting the convention 
induced by the dynamical body frame. Recall that the special xxzz gauge potentials A^^{q) 
are measurable quantities in our approach. The same holds for the angular velocity in the 
evolving dynamical body frame. 

In the static orient at ion- shape trivial SO (3) principal bundle approach the Hamiltonian 
version of the conditions 5*^' = and q^ = 0, definitory of C-horizontality and vertical- 
ity respectively, is given by Eq.( |2.32| ). Note that the definition {S^, Pfj)v = {S^, p^|g=o = 
Sq ■ Afj,{q)) of vertical rotational motion = is still valid in our dynamical body frame 
approach. By using Eg . (Pi)2|) to find ^^1^=0, we recover the rotational kinetic energy (the 
centrifugal potential) H^lf^ = ^S'^{i~^Y'^Sg of the xxzz gauge. 

On the other hand, the C-horizontal component should be determined by the condition 
Sg = 0. Since our construction requires Sg 0, we cannot utilize it. In our approach the 

measurable vibrational kinetic energy H^lf^ for S ^ non-singular N=3 configurations is 



obtained by restricting in Hrei to the value 5*^ 



of Eq. 



upon the requirement 



that the dynamical angular velocity vanishes in the xxzz gauge. From Eq. (P5|) we get 
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i f ^11 _J_ ^22 \ _|_ 



k-^{2N^-l) 

PqlPq2 



(3.32) 
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Let us remark that with this definition we get Hrei 7^ H^lf^ + Hl.lf\ differently from the 
static orientation-shape bundle result associated with the C-connection In order to get 
the theory with the Jacobi normal coordinates one has to perform our sequence of canonical 
transformations after having diagonalized kab- 



Let us end this Subsection by recalhng that in Ref. |T0[ the N=3 shape spaceQ is 
parametrized with the following configurational coordinate system {w, Wi, W2): w = p^i +Pq2 



w is the hyperradius) , Wi = p^^ — p^2, W2 = 2pgi ■ pg2 with the physical region of non- 



singular configurations defined hy ywl + W2 < w < 00, —00 < wi,W2 < 00. The variable 
w may be replaced by W3 = 2\pqi x pg2\ > {w^ > is the physical region), because 
w"^ = wf + W2 + w^. Thus, another basis for the N=3 shape space is {wi, W2, w^). A further 
basis is obtained by putting wi = wcos xcosip, Wi = wcos xsinip, W2 = wsinx, 

(x, 'ip hyperspherical angles). The quantities w and a = wl + W2 = (p^i — Pg2)^ + 4(Pgi ■ Pg2)^ 
[w2, = ^/w'^ — a] are democratic invariants, and still another basis is {w,W3,a), with a the 
angle parametrizing the democracy group SO (2). In our spin basis we have the 3 purely 

configurational shape variables pgi, Pg2 and \N\ = We can define a point canon- 

ical transformation to w = p^^ + Pg2, w^i = Pgi — Pg2? ^2 = 2pgiPg2(2iV^ — 1), and then find 
the corresponding momenta. 



C. N-Body Systems. 

Let us now consider the general case with > 4 without introducing Jacobi normal 
coordinates. Instead of coupling the centers of mass of particle clusters as it is done with 
Jacobi coordinates {center- of -mass clusters), the canonical spin bases will be obtained by 
coupling the spins of the 2-body subsystems {relative particles) Pga, Tiqa, cl = 1,..,N — 
1, defined in Eqs.( |2.7| ), in all possible ways {spin clusters from the addition of angular 
momenta). Let us stress that we can build a spin basis with a pattern of spin clusters 
completely unrelated to a possible pre-existing center- of-mass clustering. 

Let us consider the case = 4 as a prototype of the general construction. We have 
now three relative variables pgi, Pq2, Pqs and related momenta vfgi, 7fq2, vrgs- In the following 
formulas we use the convention that the subscripts a, b, c mean any permutation of 1, 2, 3. 

By using the explicit construction given in Appendix F, we define the following sequence 
of canonical transformations (we assume Sg 7^ 0; SgA ^ 0, A = a,b, c) corresponding to the 
spin clustering pattern abc 1— » {ab)c i-^ {{'^b)c) [build first the spin cluster {ab), then the spin 
cluster {{ab)c)]: 



Pqa Pqb Pqc 
T^qa T^qb ^^qc 



^^See the interpretation of the term \g^'^{q) Pfi — Sg ■ Afj,{q) Pu — Sg ■ Au{q) in Eq. (|2.34| ). 
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^It is defined for normal Jacobi coordinates pga = Sa , but it can be extended to non- Jacobi 
ones. 
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(ab)c 





ah (3b 


Olc Pc 


Pqa Pqb Pqc 


C C3 


Sqb 




T^qa TTgfc Tiqc 



Oi(ab) P{ab) l(ab) 


ttc /3c 


— * 

^(afe) Pga Pqb Pqc 




Sq(ab) Sg^ab) Sq(ab) = Sq(^ah) " ^(ab) 


C C3 


^(ab) T^qa TTgft TTgc 




«((afe)c) (3{{ab)c) 1{{ab)c) 


=? 

W({ab)c) 


=? 

7(afe) \N{ab)\ Pqa Pqb Pqc 


Sq = Sq Sg Sg = Sq ■ -/V((a6)c) 


^((a6)c) Sq(^ab) " ^(a6) ^(afe) TT^q TT^fe TTgc 



a P 'J 


— * 

W({ab)c) 


— * 

l{ab) \N(ab) 


Pqa Pqb Pqc 


Ql 52 c3 
^q ^q ^q 


^((afe)c) 


=? 7! 

^{ab) = Sq(^ab) " N(^ab) ^{ab) 


T^qa T^qb T^qc 



(3.33) 



The first non-point canonical transformation is based on tfie existence of tfie tliree unit 
vectors Ra, A = a,b, c, and of three E(3) reahzations with generators SqA, Ra and fixed 
values {R\ = l, Sa- Ra = 0) of the invariants. Use Eqs.(|37[5|), (§A^ and ( pTZl) . 

In the next canonical transformation the spins of the relative particles a and b are coupled 
to form the spin cluster (ah), leaving the relative particle c as a spectator. We use Eq.( p.l8| ) 

to define N(^ab) = |(-Ra + Rb), X{ab) = \{,Ra — Rb), S(ab) = Sqa + Sqb, Wq(ab) = Sga — Sgb- 

We get N^ab) ■ X(ab) = 0, {A^(U), ^fab)} = {Nlab).X(ab)} = { X (ab) , X (ab) } = and a new E(3) 

realization generated by S(^ab) and N(^ab), with non-fixed invariants |A''(afe)|, S(^ab) -Ni^ab) = ^(ab)- 
From Eqs. (|3.26|) it follows 



Pqa — Pqa 



\N^ab)\N^ab) + - Nlb)X{ab) 



Pqb — Pqb 
Pqc PqcRc- 



N{ab)\N(ab) 



1-N? 



{ab)X{ab) 



(3.34) 



Eq.( p.21|) defines a(^ab) and P(ab), so that Eq. ( p.22|) defines a unit vector R^ab) with S(^ab)-R{ab) = 
0) {^{ab)y ^(ab)} — 0' This uuit vcctor identifies the spin cluster {ah) in the same way as the 
unit vectors Ra = PqA identify the relative particles A. 

The next step is the coupling of the spin cluster {ah) with unit vector R{^ab) [described 
by the canonical variables S(^ab)i P(ab) 5'(^afe)] with the relative particle c with unit vector 
Rc and described by etc, Sqc, Pc, S^^: this builds the spin cluster {{ah)c). 



Again, Eq.( |3.18| ) allows to define N^(^ab)c) = \{R{ab) + Rc), X({ab)c) 



\{,R(a 



Rc) , Sq 



s 



q({ab)c) 



Sq{ab) + Sqc, 



\({ab)c) = Sq(ab) " Sqc- SiuCC WC haVC A^((a6)c) " X({ab)c) = and 
{^aaf.)c)'^f(a6)c)} = {^«afe)c) . X((af,)c) } = {X((afe)c) . X((af,)c) } = duC tO {RJ^i,),Rl} = 0, a UeW 

E(3) realization generated by Sq and N(^(^ab)c) with non-fixed invariants \N(^(^ab)c)\, Sq-N^(^ab)c) = 
Sg emerges. Eq.( p.21| ) defines a(^(ab)c) and P((ab)c), so that Eq.( p.22| ) allows to identify a final 



unit vector R((ab)c) with Sg ■ R((ab)c) 



and {Rl, 



((a6)c)' -'^({afe)c)- 



0. 



In conclusion, when Sg ^ 0, we find both a spin frame Sg, R{(ab)c), R{{ab)c) x Sg and a 
dynamical body frame X{(ab)c), A^((afe)c) x X((ab)c), A^((a6)c), like in the 3-body case. There is an 
important difference, however: the orthonormal vectors A''((afe)c) and X({afe)c) depend on the 
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momenta of the relative particles a and h through R[ab) , so that our results do not share any 
relation with the N=4 non-trivial SO (3) principal bundle of the orientation-shape bundle 
approach. 

The final 6 dynamical shape variables are = {\N(^(^ab)c)\,'l{ab), \N[ab)\, Pga, Pgt, Pqc}- While 
the last four depend only on the original relative coordinates PgA, ^ = a,b,c, the first two 
depend also on the original momenta 7?^^: therefore they are generalized shape variables. By 
using Appendix D, we obtain 

pIa = 7^^^(a,/3,7)p^A(g^P/.,^;), A = a,b,c. (3.35) 

This means that for N=4 the dynamical body frame components p^^ depend also on the 
dynamical shape momenta and on the dynamical body frame components of the spin. It is 
clear that this result stands completely outside the orient at ion- shape bundle approach. 

As shown in Appendix F, starting from the Hamiltonian H^ei{{ab)c) expressed in the final 
variables, we can define a rotational Hamiltonian -f^^e°(*(a6)c) (^^^ ^'^ ~ ^i Eqs.( [b'18| )) and 
a vibrational Hamiltonian H^^^,^^.. (vanishing of the physical dynamical angular velocity 



{{ab)c) 



0, see Eqs.( |F21]) ), but Hrei{{ab)c) fails to be the sum of these two Hamiltonians 
showing once again the non separability of rotations and vibrations. Let us stress that 
in the rotational Hamiltonian ( [I^'ISD we find tan inertia-like tensor depending only on the 
dynamical shape variables. A similar result, however, does not hold for the spin-angular 
velocity relation ( [biy] ) due to the presence of the shape momenta ^(^ab) and ^({ab)c)- 

The price to be paid for the existence of 3 global dynamical body frames for N=4 is a 
more complicated form of the Hamiltonian kinetic energy. On the other hand, dynamical 
vibrations and dynamical angular velocity are measurable quantities in each dynamical body 
frame. 

For N=5 we can repeat the previous construction either with the sequence of spin clus- 
terings abed {ab)cd ^— > {{ab)c)d) ^— {{{ab)c)d) or with the sequence abed i-^ {ab){cd) \—y 
{{ab){cd)) [a, b, c, d any permutation of 1,2,3,4] as said in the Introduction. Each spin cluster 
(...) will be identified by the unit vector R(...), axis of the spin frame of the cluster. All the 
final dynamical body frames built with this construction will have their axes depending on 
both the original configurations and momenta. 

This construction is trivially generalized to any N: we have only to classify all the possible 
spin clustering patterns. 

Therefore, for > 4 our sequence of canonical and non-canonical transformations leads 
to the following result, to be compared with Eq. ( |3.28D of the 3-body case 



PgA 



T^qA 



non can. 



a P ■y 


q^{PqA,T^qA) 


61 62 63 
^(7 


Pt^ipqA^TTqA) 



(3.36) 



This state of affairs suggests that for A^ > 4 and with Sg ^ 0, SqA 7^ 0, A = a, 6, c, 
namely when the standard (3N-3)-dimensional orient at ion- shape bundle is not trivial, the 
original (6N-6)-dimensional relative phase space admits the definition of as many dynamical 
body frames as spin canonical bases Q which are globally defined (apart isolated coordinate 



^^To each one of them corresponds a different Hamiltonian right S0(3) action on the relative 
phase space. 
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singularities) for the non-singular N-body configurations with Sqy^O (and with nonzero spin 
for each spin subclustcr). 

These dynamical body frames do not correspond to local cross sections of the static non- 
trivial orientation-shape SO (3) principal bundle and the spin canonical bases do not coincide 
with the canonical bases associated with the static theory. Therefore, we do not find gauge 
potentials as well as all the other quantities evaluated in Appendix E for N=3. 
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IV. THE CASE OF INTERACTING PARTICLES. 



Although we have defined the kinematics for free particles, the introduction of potentials 
is straightforward. 

In the case N=2, the general form of a 2-body potential is V{p^, Pq ■ 7^q, tt^) = 

V{plpq^g,^^q + 

In the case N=3, suppose there is a potential V{f]ij ■ fjhk) with fjij = ffi — ffj . Then, since 
Eqs. (|]8D and imply 

= E(ra^ - ra,)Pqa = ^ E(r.. - Ta,)Pqa[N + 

1 ^'^ 

1 2 ^ 
a=l 

+ ^[(Fi, - Tij)iT2h - T2k) + {T2i - T2j){Tlh - Tik)PglPq2i2N^ - 1), (4.1) 

it turns out that the potential is a function of p^^, p^g; (2A^^ — l)pqiPq2- On the other hand, 
if the potential is V{f]ij ■ fjhk, Hi ■ ffhk), Eqs.( |3.2(^ ) and ([C^ ) give the momentum dependence 

1,2 
a,b 



fJhk ^ i2^a^{nh - Ttk)Pqb{f^qa[N' + {-T+'+^l - N^)] + 



+ 



+ {-r^]{Sl + i-r'^^l-m). (4.2) 



By means of the same equations it is possible to study a dependence upon ■ kj. 

For > 4, even a potential of the form V{f]ij ■ ffhk)) with ffij = ffi — ffj turns out to be 
shape momentum dependent, since we have 

^ 1,..,N-1 

V {ffij - ffhk) = V[— i'^ai-'^aj)(rbh-'rbk)pqa- Pqb]- (4.3) 

a,b 

In particular, for N=4, due to Eq. ([F7|) , in the pattern {{ab)c) we get V = 

V{{ab)c)[pqa, Pqb, Pqc, \N(^(ab)c)\ , Kab) , \N(ab)\] ^{{ab)c) , ^{ab)] Sg]. 

For more general potentials V{ffij ■ ffhk, f^i ■ ffhk-, i^i ■ i^j), like the non-relativistic limit of 
the relativistic Darwin potential of Ref. [|13|, more complicated expressions are obtained. 

Finally, let us remark that both in the free and the interacting cases it is easier to 
solve the original Hamilton equations, rather than the equations for the orientation-shape 
coordinates, and then use the canonical transformation to find the time evolution of the 
dynamical body frame and related variables. 
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V. CONCLUSIONS. 



While it is possible to separate out absolute translations and relative motions in a global 
way, a global separation of absolute rotations from the associated notion of vibrations is not 
feasible because of the non-Abelian nature of rotations. 

This fact gave origin to the orientation-shape SO (3) principal bundle approach described 
in Ref. 0. While the treatment of the non-relativistic ideal rigid body is based on the 
notion of body frame 0, this and related concepts like angular velocities turn out to be 
non-measurable gauge quantities in this approach. 

In this paper we have shown that, by using geometric and group-theoretical methods and 
by employing non-point canonical transformations, it is possible to introduce both a spin 
frame and a N-dependent discrete number of evolving dynamical body frames for the N-body 
problem, with associated dynamical orientation and shape variables. Since all these concepts 
are introduced by means of a special class of canonical transformations defining the canonical 
spin bases, they, and the related angular velocities, become measurable quantities defined in 
terms of the particle relative coordinates and momenta. A notion of dynamical vibrations is 
associated to each spin basis. The spin bases are defined by a method based on the coupling 
of spins, which does not require the use of Jacobi normal relative coordinates. However, the 
Hamiltonian is not the sum of the dynamical rotational and vibrational energies: this shows 
once again that rotations and vibrations cannot be separated in a unique global way. 

Even if, in differential geometry, Darboux canonical bases have no intrinsic meaning, the 
physical requirement of the existence of the extended Galilei group with the spin invariant 
S selects a preferred class of canonical coordinates adapted to the rotation subgroup SO (3). 
Moreover, various E(3) realizations emerge whose invariants have a physical meaning. The 
adaptation to these groups selects the final class of physically preferred Darboux canonical 
bases. 

Our definition of dynamical body frames is based on the existence of Hamiltonian non- 
symmetry right actions of S0(3). We begin assuming that, like for rigid bodies, the body 
frame axes in the case N=3 be functions only of the particle relative coordinates. 

For N=3, the orientation-shape SO (3) principal bundle is trivial and we recover its 
results in a xxzz gauge. Interestingly we have found, however, that already for N=4 the it 
dynamical body frame axes depend also on the momenta so that any connection with the 
orientation-shape bundle approach is lost. 

It is an open problem whether the use of more general body frames for N=3 and > 4 
obtained by using the freedom of making arbitrary configurations dependent rotations (see 
footnote 9) may be used to simplify the free Hamiltonian and/or some type of interaction. 

It is hoped that our results may be instrumental for nuclear, atomic and molecular 
physics, since a description based on spin clustering rather than on the standard Jacobi 
center-of-mass clusterings was lacking till now. 

Let us observe that the extension of the dynamical body frames to continuous deformable 
bodies (see Ref. for an initial study of the relativistic configurations of a Klein-Gordon 
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This schematization does not exist in special relativity: the only relativistic concept of rigidity 
are Born's rigid motions |2£]. 
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field from this point of view) is a lacking piece of kinematical information and will be studied 
elsewhere. 

The fact that we use non-point canonical transformations will make the quantization 
more difficult than in the orientation-shape bundle approach, where a separation of rota- 
tions from vibrations in the Schroedinger equation is reviewed in Ref. [^J. The quantiza- 
tions of the original canonical relative variables and of the canonical spin bases will give 
equivalent quantum theories only if the non-point canonical transformations are unitarily 
implementable. These problems are completely unexplored. 

In a future paper |^ we shall study the relativistic N-body problem, where the defi- 
nition and separation of the center-of-mass motion are known to be a complicated issue. 
This problem has found a solution within the Wigner-covariant rest-frame instant form of 
dynamics [|11|]. It will be shown that concepts like reduced masses, Jacobi normal relative 



coordinates and tensor of inertia do not exist at the relativistic level. Yet, in the framework 
of the rest-frame instant form, both the oriantation-shape SO (3) principal bundle approach 
and the canonical spin bases can be defined just as in the non- relativistic case. 

ACKNOWLEDGMENTS: We are warmly indebted to our friend Michele Vallisneri for 
his irreplaceable help in drawing our frames. 
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APPENDIX A: THE EQUATIONS OF MOTION IN THE ORIENTATION-SHAPE 

BUNDLE APPROACH. 



As shown in Ref. we get the following Euler-Lagrange equations: 
i) for the orientation degrees of freedom 



dsl 

dt 



0, 
or 



^ruv 



DS, 



r 



-e^'^«[/-i'-5;]5;', (Al) 

where in the last lines D/Dt stands for the SO(3)-covariant derivative. 

For a rigid body one could use S*^ = I^'^^^^cj^ to eliminate either Sq or to from these 
equations, to obtain an autonomous system of equations (Euler's equations for a torque-free 
rigid body). This cannot be done for deformable bodies, because the relation between body 
frame spin and angular velocity involves the shape variables. Therefore, these equations are 
coupled to those for the shape variables. The motion of Sg takes place on the surface of the 

sphere = const, in body frame angular momentum space. In the case of rigid bodies, 
this motion usually follows closed orbits, but for deformable bodies the orbits will not close 
due to the coupling to the shape variables. The body frame angular momentum sphere is a 
two-dimensional surface, but counts for only a single degree of freedom, since it is a phase 
space, not a configuration space. 

ii) for the shape degrees of freedom |^ 



^ ]^gi^)^^[gi^l + gi^l - gf^ (A2) 

where (/'•'^^~^);^ = ^^^tt- — [^[f \ jC^)-!] is the covariant derivative of the inverse of the inertia 
tensor. 

While the geodesic equations (^) are 6A^ — 6 equations for the 6A^ — 6 shape variables 
g'^, Eqs.( |Al| ) are only one equation for one rotational degree of freedom. Therefore, these 



The Lagrangian is that of free particle in a non-Euclidean manifold with metric g^iu (q)- 
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equations are a system of equations for 6A^ — 5 degrees of freedom. The equations for the 
other two orientational degrees of freedom and their velocities are represented by four first- 
order differential equations: 

a) The restriction of the motion of to the surface of a sphere Sg = const., since ^ = 0; 

b) The three equations 

^R(9'') = R(e'')n(9'',r), (A3) 
at 

f -cDa UJ2 \ 

with = cDs —uJi , which can be regarded as three equations for the three Euler 

\ -UJ2 UJi / 

angles. Their solution requires a time ordering: R{t) = RoTe^^° [see Refs. [|,^]]. 

The Hamilton equations, whose content is equivalent to the Euler-Lagrange ones, are 
[Eqs.( ^^33|) are used] 

ha o v(R)oL dHrel v(R)a ~r 

dt''-~ 06-=^ dS. 



qu 



dHrel 



o dHrel / A ^ \ 
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APPENDIX B: ROTATIONAL KINEMATICS IN GENERAL COORDINATES 



In this Appendix we shall reformulate the results of Subsection D of Section II by using 
arbitrary non Jacobi relative coordinates pqa- The body frame components of coordinates 
and velocities are 



Pqa 
' r 

pqa 
Va = UJ X 



-r>TS / na\ -,s 

Pqa + -Q^q , UJ = K [U )UJ , 



Pqa = ^^ pqa + -^q^ 



(Bl) 



Let us remark that for the original particle positions r/j of Eqs.( p.8[) we get the result 

= L- + ^ Ea=l ^ai'pqa = + X {ff, - g,,) + ^^^^ ■ 

For the momenta we get 



N-l 



qa 



^qa 



J2 ^abP 



•qb 



K 



)^qa^ 



6=1 



N-l 
6=1 



N-l 



6=1 



Af-1 



^q^' 

N-l 



X ^ kabpqb{q)y + E ^'^bjrfq^, 



6=1 



while for the spin we get 



Sa 



s: 



qa 



N-l 

E pqa X TT, 
a=l 
N-l 

Eruv ~u ~v 
e PqaT^qa 



a=l 
I. .N-l 



6=1 



or e _ qr^ 
I. .N-l 

- E kabe'-'-'^Pl 

ab 



dq^' 



qaK 



E ^ab(pqa X b X p^f, + ^^(j^ 



ab 



dqi- 



(B2) 



(B3) 



By introducing the following body frame inertia tensor [{ab) means symmetrization] 

l..Af-l 



J-(g,m)= E Kbii:,){q), 



at 



sopQj. _ Q ^]^g rigid body result 6f]i = fj^dr = 6qnr + o5 x (ffi — q-nrj^T. 
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f or every pair (ab), 



fuv ( f~l \vt rut 



(B4) 



and the following potentials 



Ir^ dpqb dp^ 

l..N~l 



def 



(B5) 



ab 



we get the following gauge potentials 



(B6) 



For the body frame spin and angular velocity components we arrive at the same results 
as in Subsection D of Section II where Jacobi relative coordinates are used 

dpqb. 



Sq=Yl [Pia ■ Pqb^'' " Pqa " (^fqb + iPqa X -Q^T ^ 



ab 
I..N-1 



ah 



1..7V-1 



X-(g,m)a;^+ Kbr^^bMMabMi' 



ab 



(B7) 



I..N-1 



[X-\q,m)rSl-[X-\q,m)r{ E 

ab 

[X-\q,m)rSl - [X-\q.m)r a;{q,m)r = 

[X-\q,m)Y'S-^-A'^{q,m)r, 



For the Lagrangian we get 



1 1-^-1 . . 1 

LreliPqa, Pqa) = q E ^^bPqa ' Pqb = 7^ E ^^^bVa ' ^4 



1..7V-1 



ab 

1..N-1 



ab 



= 2 ^ p Vga ■ Pqb-I^- Pqal^ ' Pqb + 

ab 

+ ^-[p'^aX^ + Pq,x^)q +^-^qq 



dq^" 



I..N-1 



dpqa dpqb 



2 E kab[iu^)co-^'' + ^^UM)^''r + ^-^ 



(B8) 
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2 



+ 



-1 l..iV-l 



1..N-1 



+ g^lv{q, m)q^'q^\'^^ Lrei{u}{0, 9), q, q), 
where we have introduced the following pseudo-metric and metric on shape space 

l..Af-l I..N-1 



ah 



ah 



ah)tiuW - i(ah)W^(ah)aW^(ah)yVl). 



9{ah),xy\H) - ^(ab)iyj^(ab)MVy;^(a 

9t^A<l.^) = hi,u{q,m) - Al{q,m)i'''"{q,m)Al{q,m). 
The anholonomic momenta are 

dLrel 



SI 



doJr 
dLrel 

^q^' 
-1 



9nv{q, rn)q'' + Sq ■ ^^(g, m), 



d;" = [I-\q, m)rs:; - {Air\p. - Sq ■ AMq, m), 



g^'^iq, m)\p^ - Sq ■ m)], 



N-1 



Ka = e™^([X-^(g,m)]™5^^ - {A^r-'lp. - Sq ■ X])(g,m)) ^ kahP'qM + 



6=1 
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6=1 '-'^ 



(Sg, Pf,)v = {Sg, Sq ■ Ai,{q, m)), P^\g-=0 = Sq ■ A^, 

(Sg, Pi^)ch = {Q, Pix- Sq- Ai^{q,m)), p^\sr=„=Pi^- Sq- A^{q,m). (BU) 

The last lines show the phase space decomposition of generic momenta into vertical and 
C-horizontal components. 

Also the Hamiltonian takes the same form as with Jacobi coordinates 



1..N-1 



Hrel = ^ X! Kb^qa ■ TTqb = Sq ■ LJ + p^q^ - Uel 
ab 



- I [S:;{i-\q, m)rSl + {r{p, - Sq ■ A,){p. - Sq ■ X)) {q, 



ml 



^ o dHrel 

UJ = — - — 

o dH^el 

^ • 

OPi, 



(B12) 
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APPENDIX C: SOME FORMULAS FOR THE CASE N=3. 



In this Appendix we justify Eq.( ^.25D . 

To find the expression of vfga ■ N, TTga ■ x, Tiqa ■ iV x x let us consider the following quantities 

— * — * 

Sql = pql X Vfgi = Pqi[N X vfgi + f X vfgi], 
Sq2 = Pq2 X 7fg2 = Pq2[N X 7fg2 " X X ^q2], 

Sq = Sql + Sq2 = N X [PglTTqi + Pq2T^q2] + X X [Pgl^Tql - Pq2T^q2], 
Wq = Sql - Sq2 = N X [PqlTVqi - Pq2T^q2] + X X [PqlT^ql + Pq2T^q2], 

Sq-N= |iV|5j = |iV||x|[p,i7f,i - p,2V?,2] • iV X X = |iV||x|fc3, 
Sq-X= \x\Sl = -|iV||x|[p,i7f,i +p,27f,2] ■Nxx = -|iV||x|a3, 
Sq-Nxx=\Nxx\S'q = 

= N^\x\[pqlT^ql + Pq2T^q2] " X " \N \[p qlTX ql " Pq2T^q2] " N = 

= N'\x\a2 - x'\N\bi, 

Wq-N=-Sq- X, 
Wq-X=-Sq-N, 
Wg-N XX = iV^lxIbglTTgl - Pq2^q2] " X " X^ I ^^1 [Pgl???! + Pg2vfg2] " N = 

= N'\x\b2-x'\N\ai, 



PqlT^ql 
Pq2T^q2 
PglTTgl + Pq2T^q2 

PqlT^ql — Pq2T^q2 



pql ■ TTql = Pgl ^1 ■ N + TTqi ■ x], 
Pq2 ■ TTq2 = Pg2 [7?g2 " ^ " 7?g2 " X] , 
Pql ■ VTgl + Pq2 " T^q2 = 

[PglTTgl + Pg2Vrg2] " N + [pqlTlql - pq2T^q2] " X 
Pql ■ TTqi — Pq2 ■ 7lq2 = 

[PqlT^ql - Pq2T^q2] " N + [Pqingl + Pq2T^q2] " X 



|^I«1 + |X|&2, 

|iV|6i + |x|a2, 



(CI) 



where 



Oi = 


[PglVTgl + Pg2Vr<?2] 


■AT, 


02 = 


[Pgl7?gl + Pg27rg2] 


■ X, 


as = 


[PqlT^ql + Pg27?g2] 


■iVx X, 


hi = 


[PglVTgl - Pg27?g2] 


■iV, 


&2 = 


[PglVTgl - Pg2Vrg2] 


■ X, 


&3 = 


[PglVTgl - Pg2Vrg2] 


■Nxx. 



(C2) 



Eqs.(^) imply 
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1 + Pgl ■ Pq2 



[PglTTgl + Pq2T^q2] 



1 - ■ Pq2 



02 

as 

h 
b2 



1 - Pgl ■ Pq2 



1 + Pql ■ Pq2 



Sq-Nxx, 



Sq ■ X 

' \N\ ' 



1 + Pgl ■ Pg2 



1 - Pgl ■ Pg2 



[PglTTgl - Pq2'^q2] 
[PglTTgl + Pg27rg2] 



1 - Pgl • Pg2 



1 + Pgl ■ Pg2 



5„ ■ AT X X, 



W^, ■ iV X X, 



so that we get 



^ql-N 
^q2-N 

T^q2 ■ X 
Tlgl- N XX 
^q2-N XX 



. l+pql-pq2 , 5'g2-A^XX I - Pql ' Pq2 

"'^V — ^ — + p„ V — ^ — ^ 



~ Pql- f>q2 , S'gi ■ X X /l + P<yi ■ Pg2 

^9iv + ^/ ^ ■■ 



ai + 6i 

2Pqi 
Qi - ^1 

2Pg2 

02 + ^2 
2Pgl 

a2 - 62 

2pg2 

Q3 + ^3 
2P9l 

2pg2 2pg2"'' '\N\ 



1 +Pgl • Pg2 _ 'S'gl ■ X X / 1 - Pgl ■ Pg2 

2 V 2 ^ 



Pgl 



-7rg2i 
1 



1 - Pgl ■ Pg2 5*52 ■ iV X X / 1 + Pgl ■ Pg2 



Pg2 



■S„ 



X N. 



2p,i '\N\ \x\ 
1 - y iV 



(C3) 



(C4) 



This completes the study of the decomposition of the momenta Tfga on the dynamical 
body frame. 

Since Sg ■ N = Sg, Sg ■ x = sin ipcos = Sg, Sg ■ N x x = sin ipsin 7 = S'g , Wg ■ N x x = 



^yl — N"^, from Eg. (|C^) we get vfga in terms of the final canonical variables 



7fgl-iV=7fgl|iV| 



ng2 ■ N = ng2\N\ + 



2Pgl 

1-m 



2Pg2 



s 



Tfgi ■ X = vfgiVl -N^ + ]^ S'g+ - 

^Pgl 
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7rq2 • X 



Si 



Sq 



'^Pql ^ \N\ ^l-N2 



S. 



Sq 



(C5) 



For the spin quantities Wq, Sga we get 



Wq^[- 



|iv|vi -iv^ 

S„ / cos"^ ip 



sin xjjcos ipcos 7 + ^ y 1 ~ N'^sin ipsin ^]Sq + 

— N^cos Ipsin - 



[— ^l^=sm 7C0S V + - N'^cos j]SqX 



1 - iV2 

= W^X + W^N X X + W^N = W'qfr 
^Wq[\N\,^;Sq,a;SlP;Sl^], 



W„ = , cos 'ycos p + 



+ 



+ 



iSqf 



\Jl — N'^sin ^cos 



+ 



\Nr{Sq)^^Jl-N^ 



— ^ V 1 ~ N^sin 7^ j (sin asin (3 — ^cos acos (3) — 



cos^ — 



\N\S' 



+ V 1 — iV^^j sin ^{cos asin (3 — sin acos (3) , 



W 



COS J sin (3 + 



+ 



\^q> 



Vl-7V2) 



cos J 
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— p V 1 ~ N'^sin 7^ {sin acos (3 + -^cos asin j3) + 



+ 



\N\S' 



+ V 1 — N"^^ sin 'y (cos acos (3 — ^ sin asin j3), 



sis': 



\N\^l-N\S,y 



{Sqf 



+ ( . — yl ~ Jyi^jcos^cosa — 



1^1 \Sg)''^/l-m 

S'q^fiS^f^W [— 
W 



N'^cos asin 7^ + 



+ 



N\siJ{s,Y - {Slf 



1 - mSn 



-sin asin^ + 



{Sqf - {Slf 



sin acos 7^, 



\N\ ^' 



'qa 2 



(5, + {-T^W,){\N\ , S„ a; S', /3; Sl 7], 



ql _^(ql ( \a+l \^\ 
^qa - -^Wq " \- ) I 



r.2 
'-'qa 



l{s'q + {-r^'^i-N^o, 



^qa liS'q-i-r'"^^ 



\N\ 



S'q). 



We get also 



^ql — ^ql + /I ^2 



1 



e(i-N^) + (s'^y + ^{siy + 



SSL 



+ 2{i\ll-N-'Sl - 



Cl C3 
'^q'^q 



-t2 ~2 J \_ 

^q2 — ^q2 + T^I 



4P^2 



1 



e{i-N^) + {s'^Y + ^SsiY + 



- 2(eV'l - N^Sl 
1 _ 7V2 ^ 



lA^I-^l- ATS 



(C6) 
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'Pg2 Pql 



^PqlPq2 

+ (1-7V2)(27V2 



1 

7V2 



1-7V2 



+ 



+ 



(C7) 
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APPENDIX D: EULER ANGLES. 



Let us denote by a, j3, 7 the Euler angles chosen as orientation variables 6'^. 
Let /i = /2 = j, /a = ^ be the unit 3- vectors along the axes of the space frame and 
= X) ^2 = iV X X, 63 = N, the unit 3-vectors along the axes of a body frame. Then we 
have 



Slfr = R'^'{a,Pn)S'Jr 
is = {RTia, P, 7)/. = 7^/(a, P, 7)/. 



(Dl) 



There are two main conventions for the definition of the Euler angles a, /3, 7 
A) The y- convention (see Refs. |^ (Appendix B) and |^ 
i) perform a first rotation of an angle a around /s [f\ ^-^ e^ 



-sin afi + cos a/2, /s ^ 63 = /g]; 

ii) perform a second rotation of an angle (3 around 62 [^1 



cos afi + sin a/2, /2 ^ 62 
e[ = cos (3e[ — sin l3e^, 63 



-2 5 



sin Pci + cos Pe^] 



iii) perform a third rotation of an angle 7 around 63 [e[ ^— ei = cos'ye[ + sin'je^, e^ ^ e^ 
—sin^e'i + 005762]. In this way we get 



( * \ 




Nxx\ = 




\ N J 





7^(a,/3,7) 



h 



( cos^cos Pcosa — sin^sina cos^cos Psina + sin^cos a — cos sin [3 
sin ^cos (3cos a + cos ^ysin a —sin 7C0S jSsin a + cos ^cos a sin ^ysin (3 
sin (3cos a sin /3sin a cos (3 





V 




with 


tga = 




cos P = 




tgi = 




{N X xr 



(D2) 



Since and x are functions of pga only, see Eg . (13.181 ), we get {a, P} = {P, 7} = {7, a} 



0. 

B) The x-convention (see Refs. [Q, (in the text) and p|): the Euler angles 6, ip and 
ip are: i) 6 = P; ii) cos if = —sin a, simp = cos a; iii) costp = sinj, simp = —cosj. 

We use the y-convention. Following Ref. 0, let us introduce the canonical momenta p^, 
Pp, p^y conjugated to a, P, 7: {a,pa} = {P,Pf}} = {7,^7} = 1 [note that this Darboux chart 
does not exist globally] . Then, the results of Ref. imply 
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1 COS a ~ 

S = —sin ap5 H -px^ — cos actg ppa, 

sin (5 



~ stna , _ 

b — cos aps H -pj — sm actg ppa, 

sin (5 

= Pa-i 



9 



x,i . ^ COS 7 ^ ~ 

S = sm —^Pa + cos 'yctg (3p^, 

sm p 



S = cos'jpg + —^Pa - sm'yctgfJp^, 
sm p 

Sq — Pj, 

Pa — = —sin Pcos jS^ + sin [3sin jS^ + cos (33^, 
Pp — —sin aS]^ + cos aSg — sin ^Sg — cos ^S^, 
p^y — Sg — cos asin PS^ + sin asin PS^ + cos (38^. 
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APPENDIX E: THE GAUGE POTENTIAL IN THE XXZZ GAUGE. 



From the Hamilton equations associated with the Hamiltonian ( p.31|) the shape velocities 



are 



Pql 



Pq2 



\N\ 



OK 



rel 



d-Kgl 

2k~\^^,i + 2k-\2 \\N\{1 -N^)^ + (2iV2 - l)7f,2 - |iV|(l - N^)^ 

Pq2 Pq2 

dHrel _ 

,-1 



2k^'22^q2 + '2k \2 



62 
f_q 

pql 



\N\{1 - N^)^ + {2N^ - - \N\{1 - N^) — 

Pql 



rel 



2{1-N^){ 



k-\i 



+ 



k-' 



22 . 



^pI 



+ 2k-\2{l-N^) 
Since we must have 



4Pgl ^Hq2 
r2N^ 



)e + 2Vl-iV2( 



k-' 



11 



k^ 

4P^1 ^P'q2 



^PqlPq2 Pq2 Pql ' 



(El) 



7fgl|g=0 = Sq-Aql{q) = SgAl^iPql, Pq2, \N\), 
Trq2\q=0 = Sq " Aq2{q) = Sj^Ja (Pgl ' Pg2 , 

^\q=0 = Sq ■ A^{q) = S^qAliPqU Pq2, \N\), 

the gauge potential is 

•^ql,q2,(i^) — •^ql,q2,^il) = 0; 



(E2) 



k-,'k^.,^-(2N^-l)Hk-^^r 



4 



y~pr ~p2 ) 



i _|_ ^22 \ 



+ 



{l-N'^)(2N'^~l)kl 



2plp2 



\N\{1-N'^)k 



2\U-1 
12 



k-^k--^ 
'^11 '^22 



(2iv2 - mk^^Y 



k-^ 
^22 

P2 



k-lk-^-(2N-2-i){k-^y V Pi 

^{2N^-l)k^_ 
Pi 



"11 _|_ "■22 



( 



k22' {2N^ - l)k 



-1 
12 



P2 



Pi 



\N\{1 - N^)kT2' iKi , {2N'-l)ku' 



k^ik2^ - (2iv2 - mk-,^f 



pi 



P2 



pi 



P2 



2{2N^-l)k^2 



P1P2 



(E3) 
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The term in the Hamiltonian quadratic in the shape momenta allows to find the inverse 
metric 



11 



\ 



^22 



V 



iv|(i-iv^)fc7,^ 

P2 



Pi 



pi 4 

The following equations for the angular velocity 

dHrel 



"-22 



2piP2 



^ r o 



-1 rs 



'1 ^ 
CO = — 

?\r2 V 



2 A;-i,9 



+ 



12 



2 /A;-Si fc-i 



22 



aT2 V 



11 ^ fc:^ ^ fc-^2(2iv^-i) ^^^2 ^ 



L V 4p2^ 4p2^ Vp^2 p^^y 

= i-'''i<l)S', - Al{q)r{q){pu - S'^AUq)), 
2 



+ 



k ^22 



1 - 7V2 V 4p2^ 4p2^ 
|iV|Vl- Ar2' 4p^/ 



fc-\2 

2PglPy2 
22 



93 



(E4) 



0, 



X-i-(g) + ^'''^(g) At{q) 5,^Ut=o = P. 



(E5) 



allow to find the nonzero components of the inverse of the tensor of inertia in this canonical 
basis 



2 A;-i 



22 



+ 



12 



Ar2 V 4p^i 4p^2 2p,ip,2 



i-'''Hq) 



2 ,k-\i k-^ 



(■ 



22 



|JV|v'l-iV2^ '•''Si "ifl 



). 
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-1,22 



2N^ ^ Pi 



Pi 



P1P2 



\N\il-N')ki^ 



X 



-1,33 



(9) 



P1P2 

2 /A;^i 



■( 



11 



+ 



pi 

L 

22 



1 - L 



-1 L-1 



"11 



"^22 



Pi P2 



k-\,. k-' 



12 



(E6) 



The other Hamilton equations besides Eqs. (|El|) are [we do not give the exphcit expression 
of the last three of them] 



a 



rel 



1,003^7 , . 2 N/^ ^1 , ^ ^22^ 



+ sin 7)(- 



+ 



+ 



-2^ N2 ■ p2^ 



+ 



Pg2 



+ 



iV2 



2^, ^ 1,-1, 

+ (2A^2 - 1)W7)- 



PglPg2 



+ 



22 - 



+ 2fc-^2|iV|\/l-iV2(^-^) 

P<?2 Pgl 



2 Pgl Pg2 
COS 7 



+ 



,k-\i k-\2- 



2|iV|\/l-iV2 /^^i Ph 



, cos^ 7 



1 , cos^ 7 
2^ iV2 



+ sin 7 + 



1 



:)( 



k-\i , A;-i22. 



1 . A;-^ 



Pg2 



+ 



12 



1 - A^2 pg^pg2 



^,k-\, k-\2- 



Pg2 



+ 



+ 2fc-\2|iV|Vl-iV2(^-^' 

Pg2 Pgl ' 



2^ P^l 

cos^[{S,f-2{Slf]k~\, k-\,. 



2\N\\Jl-N^ 



Pql Pq2 



rel 



0, 



da 

^ ^3 ^^rel _ ^-^ .„„,,r/'c \2 I c3\2 



dt 



(97 



s^n^cos^[iS,y-iS■■:)%hi--l)i^ 

^ Pql 



11 



k-' 



22 • 



Pg2 



+ 



1 



k-' 



12 



- {2N' - 1)] 

W2 PglPg2 



M)^-(^|)^H7V'l-iV2[|(^ 

^ Pgl 



A:-^22 . 
P'2 ' 



+ 



+ 2k~\2\N\\/i-m{ 



,T[ql _ 
Pq2 Pql ' 



+ 



2|iV|Vl-iV2 
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o dHrel 
dPq2 ' 
■ ^ _dHrd 

d\N\ ^ ' 

We recover the three non-Abehan constants of motion 5"^, 5"^. 
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APPENDIX F: 4-BODY CASE 



Let us give the main steps of the construction of the spin basis in the case N = 4. From 
Eqs. (|3.34| ) we have 



Pqa Pqa 



\N(^ab)\N(^ab) + ^1 - A^(ab)X(a6) 



\N(^ab)\N^ab) - Vl-^^(afc)X(ab) 



Pqb — Pqb 
Pqc PqcRc- 

The definitions given after Eqs.( p.33D and Appendix D imply 

N — N(^l^ab)c) = -^{Riab) +Rc), 

X =^X((a6)c) = -^{Riab) - Rc), 

N XX = -^Riab) X Rc, 
4 



Pqc Pqc 



|A^|iV- Vl-A^^X 



(Fl) 



(F2) 



Let us remember that the dynamical shape variables of Eg. .331) are = { |A^|,7{afe), 
\N(^ab)\, Pqa, Pqb, Pqc}- Then, from Eqs. ( p.21| ) and ( p.24| ), for the subsystem (ab) we get 

Rab = N + X, Slab) ■ R(ab) = 0, 



^{ab) — ^qa ~r Oqh 



S{ab)[simlj(^ab) COS-i(^ab)X(ab) + siu (^ab) N (ab) X X(ab) = COS ^p(^ab)N(ab) 



X{ab) 



-{Ra - Rb) 



= COS '^(ab){sin^{^ab)S(ab) - COS ^{^ab)R(ab)) - sin'^(ab)R(ab) X S'(afe), 
N(^ab) = COStp(ab)S(ab) + siu ij(ab)R{ab) , 
N(ab) X X{ab) = sin^(^ab){sinilj(^ab)S{ab) " COS ll) {^ab)R{ab)) + COS 'y^ab)R{ab) X S(^ab), 



COS i)l^ab) 



(ab) 



sin^i^ab) 



{ab)- 



b) 



^(ab) = S[ab) ■ N(^ab)- 



(F3) 



We have moreover 
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Wq((ab)c) ■ N = -Sq ■ X, ^qi{ab)c) " X = -Sq ' N, 



C((a6)c) 
S(ab) 



\N\S^ 



1- Ar2 

\N\ 



cos^^{Sqr-{sir, 



Sq ■ X S^ 



-{Sq + Wq^(^ab)c)) = 



Slb)X + Slb)Nxx + Sl,)N 



ksl + Wiaab)c))X + + W;^^ab)c))N XX+l{S'q+ W^5(a6)c))^, 



Sqa + 'S'qb + 'S'^rc — S(^ab) + 'S'qrc — Sq(^(^ab)c) — 
S'qX + S'qN + S'qNXX, 



slab) 



d2 

^{ab) 



C.3 

^{ab) 



ny'^q I —'^q)^ 

n(Sl + i{{ab) 



oi^^ ^ \AT\ " ^i)' 



\N\ 



Sq = COS7((„b)c)V(5'g)2 - (^3)2, Sq = siu -f(^ab)c) \/ (Sq)^ - (S^)^, 



S{ab) 



'(ab) 



(Stab))' + iSlb))' + iSf^b)) 



\ 



{S^q+^iiab)c)\/l-NY + 



SI 



2S, 



{ab} 



{S'q + C^^ab)c)\/l-N^)N X X + {S'q - r^^ S^X + 
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R^ab) = \N\N+^1-N^x, 



\N\ ■ 



R{ab) X S(^a. 



b) 



is'^ + e((a.)c) VI - iv')(V 1 - Nm - \N\x) + 



\N\ 



(F4) 



'l-iV2 - - \N\ 

The final result is [5*^ = ^^^(5'^)^ - {S^Ycos P((ab)c), Sg = j^^JiSq)^ - {S^)'^sin P((ab)c)] 

N(ab) = C0S1p(ab)S(ab) + Sin Ip (^ab)R{ab) = 



^iab){S' 



Si 



\N\ ^1 



+ 



+ 



^(Siab))' 



+ \N\ 



\ 



(ab)- 



5, 



(afe) 



N + 



nSiab))' 



n^ab)iS',+^iiab)cWl-N^) 



(ab)- 



(ab) 



Nxx, 



nSiab))' 

X(ab) = COS^(^ab){sin'4)(ab)S(ab) " COS ll){ab)R(ab)) - StTl 'J (ab) R{ab) X S^ab) 



COS 7(afc) 



1 - 



^(ab) X - ^ 1 



SI 



\ S{ab) 



25, 



[ab) 



\N\ 



S{ab) 



25'(a6) 



1 - N'^sin-^i^ab) ,^2 



+ (cOS7(afe) 
+ 



1^1 



2S 



- VI - 



(afe) 



(afe) 



(ab) 



+ 



[Sg + ^{{ab)c)yi- N^)]x + 



'^S(ab) 

sin'Ji^ab) 



\ 



S(ab) "^S^ab) 



2S, 



(ab) 



msi- 



\N\ 



1 -iV2 



s')-^i-mis 



^^^^SI)])n XX. 
\N\ "^^^ 



Then, Eqs.(0) and imply 



(F5) 



Pga = ^"(«,/3,7)Pga[l^l,7(a6), \Niab)\, Pga] ^ , ^{abY, S, 
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P'qb = ^"(«,/5,7)Pg6[|^l,7(afe), \N [ab)l pqb\ t^{ab)\ S^], 

p;, = 7^'-^(a,/3,7)p;[|iV|,p,,]. 



(F6) 



This means that Pgaorb depend not only on the dynamical shape variables q'^ = {\N\ 



\N((ab)c)\,'y{ab), \N(^ab)\, Pqa, Pqb, Pqc} but also on some of the conjugate shape momenta = 
{^({ab)c},^(ab) = Sg(^ab) " N(^ab) , ^(ab) , T^qa, T^qb, T^qc} and the dynamical body frame components 
Sg of the spin. 

Therefore, a potential V{f]ij ■ fjhk) = V {pgA ■ pqs) = Vip^a^ P^b^ P^c^ Pqa ■ Pqb = (SiV^ - 
l)PqaPqb, Pqa ' Pqc, Pqb ' Pqc) bccomcs momcntum dependent. As a matter of fact the last two 
scalar products become 



Pqa ' Pqc PqaPqc 



\N^ab)\{ 



(ab) 



afe))2 



{\N\S, 



1 - msl] 



\ 



(ab)- 



b) 



+ 



+ Vi-^i^)( 



(ab) 



Pqb ■ Pqc — PqbPqc 



S{ab) 



cos -f(^ab) 



\ 



^^iab),A\N\S'q~\/l-N^Sl) 



5, 



{ab) 



b) 



\Niab)\{p^i\m 



1 - N^Sl) 



\ 



1 - (^r) - 

^{ab) ' 



(ab) 



S{ab) 



+ COS -i(ab) 



\ 



5, 



(FT) 



This spin basis could be useful when the effective potential has a negligible dependence 
on these two scalar products since ever the main bonds (for instance chemical bonds in 
molecular dynamics) do not involve them. 

By using Eg. (|C5|) , for the canonical momenta we get [\N\ = \N({ab)c)\, S^ab) given by 
Eq.(0)] 



qa 



TT. 



qa 



\Niab)\ 



1 - , 



2p, 



qa 



(Siab))^ - {^{ab)ysin^{ab) + ^{ab)Jl- N^ab))\N{ab) + 



TT, 



qa 



qa 



{S{ab)y - {^iab)ysin-f^ab) + ^{ab)Jl- A^Jfe) 



X{ab) 



2P 



{S{ab))^ - {^(ab)} 



-,COS -f(^ab) 



(ab) 



qa 



--)N{ab) X X(a6), 



I _ I 

^qb\N{ab)\ + ^ ^ ^"''^ (\/ ('g(afe))^ " (ab)? siu {ab) " ^{ab)^l- iVjfc)) 



(Siab)? - (fi(a6))^Sm7(„b) - ^{ab)Jl " A^Jfe) 



(ab) 



X{ab) 



2pqb 



{S(ab)y - {^(ab)} 



; COS 'J(ab) 
' \Niab)\ 



^(ai 



{ab) 



1 - ^U) 



)N{ab) X X(afe), 
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TT, 



qc 



2p, 



(5,^-e((a5)c)Vl-iV') 



qc 



^^^^ 



N - 



X - 



^qc 



1 



TT 



(jra (jra 



4p: 



2 /I 7\T2 



qa 



+ [Sin li^ah) + ^ J((>g(ab)) - {^iab)) ) + "^-^ 



+ 



(afe) 

^{ab)COS ^(^ab) 



^2 ~2 I 
T^qb = '^qb + 



(ab)j 



, / . 2 COS 7(afe) 

+ [sin 7(a6) + — 



{{S{ab) f - (^(a6))^) + 



fi{afe)COS7(ab) 



\Niab)\^l-Nl,^' 



TT 



9C "gc ' ^^2 



e(W)(i-^) + (^.) + 



qc 



Ar2 



+ 



-2(e({a6)c)Vl-iV'5, 



Cl C3 



\N\^l-N^ 



(F8) 



To evaluate Tr^c one defines an auxiliary variable T[(^ab) such that S(^ab) = R(ab) x 7f(a5) [its 
component along R(ab) is arbitrary; remember that R{ab) ■ S(^ab) = 0]. Then we use the 
equations Sq = S(^ab) + Sqc and Wg = S(^ab) — Sgc (with Sgc = Pqc x Tr^c) to extract the form 
of TTqc like in the 3-body case [Eg. (pB])]. 

Considering a set of Jacobi coordinates pqA, A = 1, 2, 3, with reduced masses we get 
the following expression for the Hamiltonian to be compared with Eq.( p.34| ) [|A^| = |A''((afe)c) |] 



rel((ab)c) 



TT^ TT^ TT"^ 

qa qb ^ qc ^ 



2/ic 2/i<j 2/i;, 2/i, 



+ TZ {sin^ -f(^ab) + 



cos 



l{ab) ■ 



{^{ab))" 



^^{ab) 



1 - 



+ 



^ r 



ab)Sin^^ab)]ll- N^ab) 
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^{ab)COS 7(afe) 



+ 



I^(a6)| 
1 



/V2 



i({afe)c)l 



qc 



Cl C3 



\2 



71 = 7l{\N\,^(^(^ab)c),^{ab); Sq) = (5'(a6)) - {^{ab)] 
_ 1 

~ 4 



+ 



(5,^+e((a6)e)Vl-iVT + 

5-1 

{^{ab)] 



g 

liV 



1- 



(F9) 



Due to the presence of the function TZ, this Hamiltonian is different from that of the 
static orientation-shape bundle approach given in Eq.( p.34] ) with N = 4: note that it is not 
a polynomial of second order in the dynamical body frame components of the spin. Yet, the 
rotational part of the Hamiltonian, determined by putting the dynamical shape velocities 
equal to zero {q^ = 0), can be shown to have the standard form of the rigid body case, with 
a non-standard dynamical inertia-like tensor depending only on the shape variables. 

To show this result, let us consider the Hamilton equations for the velocities 



d 
dt 



Pqa 
Pqb 

Pqc 
{(ab)c) I 



o_ dHrel{{ab)c) 


_ T^qa 


dllqa 


Pa ' 


o dHrel{{ab)c) 


_ TTgft 


dTTqb 


Pb ' 


o dHrel{{ab)c) 


_ TTgc 




Pc ' 


o_ dHrel{{ab)c) 




d^{iab)c) 





2 

1 , 1 



(5,^ + e((a.)c)Vl-A^')x 



8 PaP^ 



+ 



qa 



PbPlb 



){sin^ 'yi^ab) + 



COS^ Ijab 
^^{ab) 



+ 



+ 



16 Papl 



qa 



PbPlb ' 



){Uab)dl- Nf^^^sinj^ab) 



^{ab)COS ^(^ab) 



(ab) 



7^ 



qc 



l{ab) 



rel{{ab)c) 
dQ{ab) 
1. 1 

+ 



qa 



PbPqb 



2 J^iab) 
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1 - ^fa^) 



cos 



I'iab) 



-)- 



COS -'i\ah) 



(^(ab))' 



8 IXaPla |A^(a6)IVl-A^J,6) 
o_ dHrel{{ab)c) _ 

11 1 / ^» 

O l^aPqa PbPqb * ^ ^ 

If we put all the shape velocities equal to zero we get 

i) pqA = 0, A = a, 6, c, 

7rgA|<j=o = 0, /l = a, 6, c, 



7^ 



^0 ^I^(a6)|=0, 





■f i{iab)c)\q= 


( 






v/l-iV2 


1 

^ 1 


1 



- (^^(a6)|q=o)^) 

(ab) 



1/2 



1- at; 



2- Stn-f(^ab) 



-^(ab) I 



g=0) 



7(a6) = 0, 

16 /^^P^.^^^"'^""^^^-^MiV(„,)|(l-Afj,))^e(a.)|,^o 



1 I 1 



cos 



l{ab) 



1 



f^(a6)|9=0 4cOS7(a6) 



= 0, 



4 



^{ab) \q=0 = X±{q) i{ab)\q=0, foT ^(ab)\q=0 7^ 0, 



Ai6P^(,\2 S«n 7(afe) COS 7(a5) 
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cos 



l{oh) 



1 - K^) 



± 



cos" 7(ab) 



COS" 7(„fe 



and then 



dt 



( 



1 1 



+ 



COS^(ab) ^{ab)\q=0 



. 2 , COS 7(afc) 

s^r^ 7(a&) H + 



s^n 



7(a6)( 



'^PcPqc 

Therefore we get 

C((a6)c)|q=0 = gdQ) 



IJ-aPqa WPg6 



- sin 7(afe)) 



l/2n 



(F13) 



(F14) 



1 



/ . 2 , COS 7(afe) 1 / tiaPla Pbpi^ \ 

[sm 7(a6) + — l( 1 , 1 J 



12- H 3- 



s^r^7(a6) szn7(„b) 



l^(a^)IVl-^(a;>) 



X 



COS" 7(ab) 1 / MaPia MbP^6 \2 



s^r^7(„^,) 



s^r^ 7(a&) 



cos7(ab)X±(g) 



1 



|jV(ai,)l\/l-iVf„», 



MaP^a WP^ 



Then, from ii) and iii) we get for ^(ab)|g=o 7^ 



M4=o = ±^([(1 + 9MW^-N']iSlf + ( J- 



1^1 yr^' 



1/2 



X 



s^r^ 7(a6) 



^(o6)|?=0 = -'^±(?)C(o6)|q=0- 



(F15) 



(F16) 
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The non vanishing shape momenta = p^\ci=o do not have the form Sg ■ A^{q) of the 
static orientation-shape bundle approach. A discussion of the various solutions will be done 
elsewhere. Here we are mainly interested in determining the form of the rotational kinetic 
energy. From Eq. (|F16|) we get 

\2 d.f^f' 



b)\g=0) 



9 9' 



,2 ' 



^ a/1 - N, 



(ab) 



^liab)iQ)i 



9=0 



1 - Kb) 



- -V ■5^'^ 7(afe) 



^(a6)|ij=0; 



e(%.)e)l.=0(l-iV^) + ^ + (^,Y + 



ATS 



-2(e(a, 



,6)c)|g=0 



1 - A^2^2 



Cl C3 



{SI? 

l-iV2 

def 



|iV|Vl- 



Therefore, from Eq.(^^ we obtain 



(rot) 



rel 



rel I q=0 



(FIT) 



J-(g)=T[:,)(g)+Tr(g) 



1/ ^ , ^ \\( ■ 2 COS l{ab).^rs( N , 

1 2-^^(-)^-(^)^x-,„,(,)-(i-ivf.jxru(.) 



(F18) 



where is an inertia-like tensor depending only on the dynamical shape variables and 

Tl'^i^q) is defined in the last line of Eg. ( It*' 17]) . 

On the other hand, the dynamical angular velocity evaluated in the dynamical body 
frame is 



^l{ab)c) 



^\{ab)c) 



rel{{ab)c) 



dSgr 



16\N\ 
/ 1 



/I , 1 W . 2 , C0S^7{ab)\ 



+ 



^{ab)Sin^^ab)^Jl- N^ab) 



2Vn 
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4/Xep2JiV| ^ \N\ 
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+ 



\N\ 
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((a6)c) 
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2Vn 
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iV2 1^1 x/l-iV2 



+ 



(F19) 



We see therefore that for > 4 in the spin bases the spin- angular velocity relation is 
different from the result ( 2.33|) of the static orient at ion- shape bundle approach. 

If we define the dynamical vibrations as those corresponding to the vanishing of the 
measurable angular velocity (there is no need of connections like in the static orientation- 
shape bundle approach), we get 



r-l 
^q 

Si 



Si 



0, 



=0 = fl{.q)i(ab)c) + f2{q)^{ab) + f3{q)^(ab) ^= /''(g)p/., 
so that the vibrational Hamiltonian is 

a,b,c ~2 
_i 

A 



(F20) 
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COS^ l{ab) 
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(F21) 



Let us stress that the Hamiltonian is not the sum of the rotational and vibrational parts: 

[(rot) I Trimb) 

^rel{{ab)c) -'-'reZ((a6)c) • 



Hrel{{ab)c) 7^ 
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By calling 9" the Euler angles a, f3, 7, the other Hamilton equations become 
is the inverse of A(^)(")^ after Eq.(^^); see Appendix A for the comparison with the static 
orient at ion- shape bundle approach] 

na o ^(R)a dHrel{{ab)c) iy(R)a ~,r 

= X' ' ^f3=r) -A^ {I3=r)^{iab)c), 

d ^ 9Hrel{{ab)c) sru n dHrel{{ab)c) ^ 



■'qu 



rsu Q ^ u 

^qs^{{ab)c)^ 



d _ o dHrel({ab)c) . , 

-i:^qA = o , A = a,b,c, 

at dpqA 

d ^ dHrel{{ab)c) 

^^"'^'^'^^ " oinT' 

d ^ dHrel{(ab)c) 

Uiab) = (F22) 

dt d'J^ab) 

For N=4 Ref. [|1^] gives the static shape coordinate basis w = p'^^ + + Pqs 

the hyperradius), Wi = ^{pgi - Pg2)? ^2 = V^pqi ■ Pq2, W3 = \pipq2 ■ pq?., ^4 = \Pipqz ■ Pql, 
= \ [— — + ^Pqs] ]V = Pqi' {pq2 X Pq^) must be added to solve the problems of shapes 
connected by parity. By using Eqs.(p^) these variables can be reexpressed in our spin basis, 
in which they depend on the dynamical shape variables, on the dynamical shape momenta 
and also on the dynamical body frame components of the spin. Therefore, to recover locally 
the above N = 4 canonical basis of the static orientation-shape bundle approach, a non-point 
canonical transformation on all the variables is needed. 



Note that there is in addition wq = \/^\J\pqi x py2P + \Pq2 x Pgsp + IPgs x Pqip > 0, 
with w"^ = J2f=iwf. There are 3 democratic invariants w, a = J2^=iwf, and b [see (4.25) 
of the paper]. Due to the absence of the degeneracy submanifolds of planar configurations, 
collinear configurations and 4-body scattering point, the physical region has a complicated 
definition based on 3 eigenvales Ai > A2 > A3 > and is topologically B = R x (S^ — P), 
with P RP^ being the 2-dimensional projective plane. 
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